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ANALYSIS OF FATIGUE, FATIGUE -CRACK 
PROPAGATION, AND FRACTURE DATA 

By Carl E. Jaske, Charles E. Feddersen, 

Kent B. Davies, and Richard C. Rice 

Battelle*s Columbus Laboratories 
SUMMARY 

Analytical methods have been developed for consolidation of fatigue, fa- 
tigue-crack propagation, and fracture data for use in design of metallic aero- 
space structural components. To evaluate these methods, a comprehensive file of 
data on 2024 and 7075 aluminums, Ti-6A1-4V, and 300M and D6AC steels was estab- 
lished. Data were obtained from both published literature and unpublished 
reports furnished by aerospace companies. Fatigue and fatigue-crack-propagation 
analyses were restricted to information obtained from constant-amplitude load or 
strain cycling of specimens in air at room temperature. Fracture toughness data 
were from tests of center-cracked tension panels, part-through crack specimens, 
and compact-tension specimens. 

Both fatigue and fatigue -crack-propagation data were analyzed on a statis- 
tical basis using a least-squares regression approach. An arc-hyperbolic tan- 
gent function was used to relate the independent variable to the dependent vari- 
able. For fatigue, an equivalent strain parameter was used to account for stress 
ratio effects and was treated as the independent variable, and cyclic fatigue 
life was considered to be the dependent variable. An effective stress -intensity 
factor was used to account for the effect of load ratio on fatigue-crack propa- 
gation and was treated as the independent variable. In this latter case, crack- 
growth rate was considered to be the dependent variable. 

Smooth-specimen and notched-specimen fatigue data were treated separately. 
Notched-specimen results were analyzed using a local stress-strain approach to 
account for fatigue damage at the notch root. Data for various types of notches 
and theoretical stress -concentration factors were consolidated by using a com- 
puted fatigue-strength reduction factor. Both the cyclic and monotonic stress- 
strain curves were employed in calculating the local stress-strain response from 
nominal loading information. 

After computing mean fatigue and crack-growth curves by least-squares 
regression, tolerance -level curves were determined. Lower-level tolerance curves 
for 90 and 99 percent probability of survival with 95 percent level of confidence 
were determined for each fatigue curve. Two-sided tolerance bands for 90 and 99 


percent probability with 95 percent confidence were determined for each mean 
crack-growth curve. 

Fracture toughness data were tabulated for a particular material and speci- 
men thickness in terms of average values at various temperatures and panel widths. 
Apparent, critical, and onset fracture toughness indexes were used in this tabu- 
lation. 


INTRODUCTION 

Recent experience with modern aerospace structures has emphasized the impor- 
tance of considering both fatigue and fracture in the design and service per- 
formance of aircraft. A structural member may fracture at loads well below the 
nominal yield strength of the material if it contains a critical-size flaw. In 
some instances, such flaws may be introduced into the structural material by 
manufacturing processes. However, in most cases, flaws will become critical by 
growing from smaller flaws or from unflawed areas of stress concentration. This 
type of growth occurs by fatigue processes during cyclic loading of the structure. 

Reaching total fracture under cyclic loading involves both fatigue-crack 
initiation and propagation. As shown in figure 1, crack-initiation life can vary 
considerably, depending upon the definition of a crack. The wide range of sizes 
(10 ^ to 10”*^ inch) considered to be cracks by various investigators causes an 
ill-defined area of overlap between initiation and propagation. in most fatigue 
tests of small specimens of virgin material, the initiation phase is generally 
considered to be a more significant portion of cyclic life than the propagation 
phase. Fatigue -crack-propagation data are usually obtained from precracked or 
flawed specimens. 

In this program, fatigue data from uncracked smooth or notched specimens 
were treated separately from fatigue -crack -propagation data from precracked speci- 
mens. It was assumed that the total number of cycles to failure normally reported 
in fatigue tests of simple specimens was a reasonable approximation of the number 
of cycles required to initiate an engineering size flaw. Crack-propagation infor- 
mation was obtained from studies where cyclic crack growth was measured using a 
precracked sample. Fracture at a critical load level or flaw size also was 
treated separately. 
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For conventional static properties of metals and alloys, extensive design 
allowables information is available in documents such as MIL-HDBK-5B (ref. 1). 

For fatigue, fatigue-crack propagation, and fracture data, however, design 
allowable values are usually not available and the data are presented in terms 
of typical or average values. 

Part of the problem for fatigue and fatigue-crack propagation is that these 
behaviors are influenced by a wide range of parameters that include cyclic stress 
mean stress, cyclic frequency, temperature, environment, product form and orien- 
tation with respect to loading, structural geometry (size, shape, and notch con- 
figuration), metallurgical and surface effects associated with heat treatment, 
microstructure, and machining practices. Most aerospace companies tend to gen- 
erate data for a limited number of these many variables to fulfill specific local 
design needs. Much of this information is retained within each company, and that 
which becomes available in open literature is often digested in accordance with 
particular theoretical considerations and analytical procedures endemic to a 
given organization. Since these considerations and procedures vary among com- 
panies, it is difficult to affect a systematic consolidation of such data. As- 
sessment of fatigue and fatigue -crack-propagation data is further complicated by 
the fact that there have been no standard methods for these types of testing. 
Recommended procedures for fatigue testing have been published recently (ref. 2). 

Many of the aforementioned problems also influence fracture results. Stan- 
dards for obtaining plane-strain fracture-toughness information have been devel- 
oped recently (ref. 3). However, major differences in testing and analysis pro- 
cedures still exist for plane-stress and transition-thickness conditions. 

As a result of extensive visits and discussions with all major aerospace 
companies, it became evident to personnel at Battelle's Columbus Laboratories 
(BCL) who are responsible for maintaining MIL-HDBK-5B (ref. 1) that the major 
deficiencies in the Handbook were in the important areas of fatigue, fatigue - 
crack propagation, and fracture. This realization led to the initiation of a 
research program at BCL under the sponsorship of Langley Research Center. Work 
was directed toward systematizing and consolidating available fatigue, fatigue - 
crack propagation, and fracture information on 2024 and 7075 aluminum alloys, Ti- 
6A1-4V alloy, and 300M and D6AC steels. It was considered imperative that the 
analytical procedures be compatible with statistical methods of data presentation 
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Similar approaches were used for both fatigue and fatigue-crack propagation, 
as illustrated in figure 2 . The logarithm of fatigue life was the dependent 
variable in both cases. An equivalent strain parameter similar to that suggested 
by Walker (ref. 4) and Smith, et al. (ref. 5) was used to account for stress 
ratio effects and was treated as an independent variable in the fatigue analysis. 

A similar effective stress-intensity factor (ref. 4) was used to account for 
stress ratio effects and was treated as the independent variable in the fatigue- 
crack-propagation analysis. 

Fatigue-crack propagation is more complicated than fatigue because differ- 
ent life curves (fig. 2) are obtained for each different state of initial damage. 
Thus, fatigue-crack-propagation results are usually presented in terms of crack- 
growth rate as shown schematically in figure 3. The layering of rate data as a 
function of stress ratio can be accounted for using the effective stress-intensity 
concept mentioned above. 

Treatment of fracture data was limited to tabulation and graphical summary 
of information in terms of three indexes of fracture toughness. 

Primary emphasis of the fatigue work was on data for 2024 and 7075 aluminum 
alloys, with a secondary effort directed toward annealed Ti-6A1-4V alloy and a 
high-strength steel. The fatigue-crack propagation and fracture work was limited 
to data for the 2024 and 7075 alloys and 300M and D6AC steels. 
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SYMBOLS 


A. 

1 


a 

da/dN 

Aa/AN 

a. 

1 

B 

C 

Cl ,Cj2 


D 


DD . , . . 
J(i) 


e 

Ae 


e 

a 

e,0),%(2> 

e 

max 

Ae 

P 

K 

AK 


K 

app 

K 

c 


K 


eff 


mean stress coefficient 
ith regression coefficient 
crack length, mm (in.) 
fatigue crack growth rate 
approximate crack growth rate 
ith value of crack length, mm (in.) 
stress amplitude coefficient 
Paris * coefficient 

regression coefficient in arc- hyperbolic tangent relation 
crack half length, mm (in.) 

crack half length at unstable fracture, mm (in.) 

initial crack half length, mm (in. ) 

multiplicative regression coefficient 

the jth order divided difference at point i 

elastic mo du lus , MN/m^ (ks i ) 

nominal strain 

nominal strain range 

nominal strain amplitude 

specific values of e used to define stress-strain curve 

a 

maximum nominal strain 

plastic nominal strain range 

stress intensity factor, MN/m^^^ (ksi/in. ) 

range of stress intensity factor, MN/m^/^ (ksi/TrT. ) 

apparent fracture toughness, MN/m^/^ (ksi/irT. ) 

critical fracture toughness, MN/m^/^ (ksi/in. ) 

effective stress intensity factor, MN/m^/^ (ksi/in. ) 
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K, 


Ic 

K. 


K 

max 

Ko 

K 

q 

K 

e 

K 

a 

Ki,K3 

k 

L 

M 

m' 


m 

N 


N, 


N. 


n 


n. ,n. 


P 

max 


Q 


q 

R 

R2 


R 


m 


2 


fatigue strength reduction factor 

mode I critical fracture toughness, MN/m^ (ksi/Tn. ) 
ith value of stress-intensity factor, MN/m^ (ksi/TiT. ) 
maximum stress-intensity factor, MN/m^ (ksi/THT) 
threshold stress-intensity factor, MN/m^ (ksi/in. ) 
toughness at onset of crack extension, MN/m^ (ksi/in.) 
theoretical stress concentration factor 
strain concentration factor 
stress concentration factor 
strength coefficients, MN/m^ (ksi) 

Stulen coefficient 
regression coefficient 
Elber optimization coefficient 
value of slope of function $ 

Walker exponent 

number of cycles 

number of cycles to failure 

ith value of number of cycles 

sample population 

Paris* exponent 

strain hardening exponents 

maximum load, KN (kip) 

5 percent secant offset load , KN (kip) 
plasticity corrected shape factor 
degree of polynomial 
stress ratio 

proportion of variation explained by regression equation 
modification of R^ to account for degrees of freedom 
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S 

AS 

SSD 

S 

a 

S^(l),S^(2) 

s 

c 

S 

eq 

S 

m 

S 

max 

S . 

min 

S 

n 

S 

o 

S 

op 

s 

T 

TSS 

TUS 

TYS 

t 

u 

V 
w 
X 

X 

Y 

Yi 
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notch root radius, mm (in.) 
nominal stress, MN/m^ (ksi) 
nominal stress range, MN/m^ (ksi) 
residual sum of squares 
nominal stress amplitude, MN/m^ (ksi) 

specific values of S used to define stress -strain curve 

a 

maximum nominal fracture stress, MN/m^ (ksi) 
equivalent nominal stress, MN/m^ (ksi) 
mean nominal stress, MN/ra^ (ksi) 
maximum nominal stress, MN/m^ (ksi) 
minimum nominal stress, MN/m^ (ksi) 
net stress or uncracked section 

nominal stress at onset of crack extension, MN/m^ (ksi) 

crack closure nominal stress, MN/m^ (ksi), S >0 

’op 

standard error of estimate 
thickness, mm (in.) 

regression (or total) sum of squares 
tensile ultimate strength, MN/m~ (ksi) 
tensile yield strength, MN/m^ (ksi) 

Students* t multiplier 
confidence level 
variance 

panel width, mm (in,) 

independent variable - tanh"^ tanh“^ [§(K^ff)] 

ith value of independent variable 
mean value of X 

dependent variable - log Nf or log da/dN 
ith value of dependent variable 
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P 

e 

Ae 
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e 
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eq 

e 

max 

Ae 

P 

e 
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P 

a 

Aa 


m 


max 


cj . 

mxn 


mean value of Y 

tolerance limit at a level of confidence (u) for a given number 
of degrees of freedom (v) 

mean stress exponent 

stress amplitude exponent 

local strain 

local strain range 
translated local strain 
local strain amplitude 

lower limit of inverse hyperbolic tangent function 
equivalent local strain 
maximum local strain 
plastic local strain range 

upper limit of inverse hyperbolic tangent function 

degrees of freedom 

notch analysis material constant 

local stress, MN/m® (ksi) 

local stress range, MN/m^ (ksi) 

translated local stress, MN/m^ (ksi) 

local stress amplitude, MN/m^ (ksi) 

mean local stress, MN/m^ (ksi) 

maximum local stress, MN/m^ (ksi) 

minimum local stress, MN/m^ (ksi) 

value of intercept of function f 
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ACQUISITION, COMPILATION, STORAGE, 
AND RETRIEVAL OF DATA 


To implement the evaluation of existing fatigue, fatigue-crack-propagation, 
and fracture data, it was necessary to make an extensive survey of the literature 
and of aerospace companies that might have unpublished internal reports. A 
computerized system was developed to compile and store data obtained from this 
survey. Data from more than 120 reports and documents were acquired, compiled, 
and stored. 


Data Acquisition 

Information was taken both from the open literature and from company reports. 
Applicable reports were obtained from the technical files of the Metals and 
Ceramics Information Center (MCIC) located at BCL. Throughout the program, new 
reports, acquired by MCIC, were screened and added to the data base when applica- 
ble. In order to obtain as much recent information as possible, additional liter- 
ature searches were obtained from the National Aeronautics and Space Administra- 
tion (NASA) (refs, 6 and 7) and the Defense Documentation Center (DDC) (refs. 8 
and 9), In addition, pertinent reports obtained through the MIL-HDBK-5 (ref. 1) 
program were used. 

Internal reports from aerospace companies and unpublished data were obtained 
from various laboratories that conduct fatigue, fatigue-crack-propagation, and 
fracture research, A letter was prepared and sent out to 89 selected members of 
the American Society for Testing and Materials (ASTM) Committee E09 on Fatigue, 

A similar letter was also sent to 18 members of the ASTM Committee E24 on Frac- 
ture Testing of Metals and to 46 members of the MIL-HDBK-5 Coordination Group. 
Positive responses were received from about 40 percent of those surveyed. 

Pertinent data from the responses were entered into the data storage files. The 
type of information that was requested in these letters is summarized in the 
following three sections. 

Basic Fatigue Information .— For the alloys of interest (2024 and 7075 
aluminum, Ti-6A1-4V, and high-strength steels), fatigue data were desired from 
axial-load tests of simple specimens that reflect basic material behavior. This 
requirement excluded joints or components but included both notched and unnotched 
data, where notch configuration and severity were variables. Data for cyclic 
12 


lives ranging from >10® to <10”^ cycles, both strain and load-controlled test data, 
and variable stress ratio (or mean stress) data were of interest. 

Basic test data were desired; i.e., tables of stress or strain versus life- 
time. In cases where crack initiation was determined, this information (and the 
initiation criterion employed) was desired. For tests involving cyclic plasticity, 
cyclic stress-strain information in the form of stress and strain as a function of 
loading history were needed. 

In addition to the fatigue data, correlative information concerning specimen 
geometry and fabrication, material product form, dimensions and processing, test 
techniques and controls, laboratory environment and mechanical properties were 
also desired. The latter information was required to aid in making decisions 
about pooling various samples of data. 

Fatigue -Crack -Propagation Information . “Fatigue-crack-propagation data were 
desired for center-cracked panels (in a variety of widths), part-through-cracked 
or surface -flawed specimens, compact- tension specimens, and double-cantilever- 
beam specimens. It was useful to have data that delineated crack initiation 
cycles from a geometrically known starter flaw as well as initial propagation 
data from it, if such information was available. Delineation of the stress cycle 
employed for each test, as well as test frequency, was necessary. In some cases, 
multiple tests were conducted on a single specimen such that propagation occurred 
on successive crack-growth segments under different cyclic-stress conditions. 

Each of these conditions was considered as a single test in the analysis, and 
the conditions needed to be described. 

Basic test data again were desired; i.e., tabular displays of crack size 
versus cycles. For each specimen, the associated test stress cycle description 
was given. 

In addition to the basic crack-propagation data, correlative information as 
described for fatigue data were necessary. 

Fracture Information . — The fracture data collection was more complex in that 
there was a thickness dependence on fracture toughness that was of greater signif- 
icance than for fatigue and fatigue-crack propagation. This thickness dependence 
influenced the mode of fracture (such designations as slant, transition, and flat 
fracture descriptions were used) corresponding to plane stress, transitional 
stress, and plane-strain fracture toughness. A variety of tests have been 
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employed, only a few of which are standardized (ref. 3), Thus, test specimen 
description and test techniques had to be delineated carefully and completely. 
Specific tests for which data were desired included the center-cracked panel, 
part-through or surface -flawed specimens, compact -tens ion specimens, double- 
cantilever-beam specimens, and notched-bend specimens, 

Basic test data were needed rather than fracture toughness values. These 
included original crack length, critical crack length, ultimate load or stress, 
and load or stress at which slow stable crack growth initiated (presented in 
tabular form). Load -compliance records were obtained when possible. 

As stated before, correlative information was desired with particular 
emphasis on test methods and techniques. 

Data Recording and Storage 

Information used in this program was stored in a format for computerized 
analysis. Detailed data has been recorded on punched cards for use at BCL. 

These data were transferred to magnetic tapes and forwarded to NASA Langley Re- 
search Center, To help document the encoded data, a short abstract was prepared 
for each report from which information was taken. This abstract summarized 
briefly the type of data encoded along with correlative information not recorded 
in the data file. The check list shown in figure 4 was used in preparing these 
abstracts. Sample abstracts are presented in figures 5 and 6. A complete set 
of abstracts was sent to NASA Langley along with the magnetic tapes. Each source 
from which data were taken and an abstract prepared are listed in Appendix A. 

Each source was assigned a unique reference number when it was added to the data 
base . 

The basic medium for recording the fatigue, fatigue -crack-propagation, and 
fracture data collected and compiled on this program is the standard 80 column 
computer punch card. Data card file sequences and formats which have been se- 
lected to provide a consistent procedure for encoding these data are described 
in the following subsections. 

Each data file may contain up to four basic types of cards depending on the 
type of information being recorded. These card forms are 

Card 1: Title or lead card, identifying test and material 

Card 2: Subtitle card, containing supplementary testing, composi- 

tional, or processing information where desirable 
Card 3: Data card describing specific test parameters and results 

Card 4: Crack growth card listing cycle count and crack size. 
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DATA SOURCE ABSTRACT CHECKLIST 


For Each Report From Which Data is Obtained 
Check for, and Record, the Following Items. 

General Report Information 

(1) Reference Number. 

(2) Materials. 

(3) Authors, Title, Publisher/Source, Publication Date. 

Test Information 

(1) Type of Test (Fatigue, Fatigue-Crack Propagation, Fracture), Summary 
of Report Abstract. 

(2) Type of Test Machines, Load or Strain Control? 

(3) Number of Specimens. 

(4) Stress Ratios. 

(5) Test Temperature and Environment. 

(6) Test Frequencies. 

(7) If Fatigue-Crack Propagation — 

(a) Plane Strain or Plane Stress? 

(b) Basic Data or "Digested" Data? 

Specimen Data 

(1) Melting Practice/Heat Treatment of Specimens. 

(2) Ductility. 

(3) Fabrication Methods. 

(4) Surface Finish. 

(5) Specimen Dimensions. 

(6) Chemical Composition. 

(7) Tensile Properties (TYS, TUS, Reduction of Area, Elongation, Elastic 
Modulus). 

(8) Are There Stress-Strain Curves or Data? Are They Monotonic or 
Cyclic? 

Figure 4. — Checklist used for preparation of report abstracts. 
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REFERENCE NUMBER 3 

Materials: 2024-T3, 7075-T6, 4130 Steel 

Grover, H. J., Bishop, S. M., Jackson, L. R., "Axial-Load Fatigue Tests on Notched 
Sheet Specimens of 24S-T3 and 75S-T6 Aluminum Alloys and of SAE 4130 Steel with 
Stress-Concentration Factor of 5.0", NACA TN 2390, Battelle Memorial Institute, 
June (1951). 


Test Information 


(1) Fatigue Tests i Axial-load fatigue tests were conducted on notched specimens 
of three sheet materials with one stress concentration factor and four mean 
stress levels. 

(2) Type of Test Machine: Krouse direct repeated-stress test machine. 

(3) Number of Specimens: 49/2024-T3, 47/7075-T6, and 42/4130. 

(4) Stress Ratio: R = -1.0 to 0.70. 

(5) Test Temperature and Environment: Tests were conducted at room temp>erature 

in air. 

(6) Test Frequency: 1100 - 1500 cpm. 


Specimen Data 

Cl) Melting Practice/Heat Treatment: Not specified. 

(2) Ductility: Not specified. 

(3) Fabrication Methods; Specimens were machined from 0.09 inch thick 2024-T3 and 
7075-T6 aluminum and from 0.075 inch thick 4130 steel. Notches were cut in a 
series of machining cuts. 

(4) Surface Finish: Specimens were electropolished . 

(5) Spe cimen Dimensions: Gross length = 15.5 inches, gross width =2.25 inches, 

net width = 1.5 inches, root radius = 0.03125 inch (K^ = 5). 

(6) Chemical Composition: Not specified. 

(7) Tensile Properties : 



TYS, 

TUS, 

Elong. , 

Material 

ksi 

ksi 

7o 

2024-T3 

54.0 

73.0 

18.2 

7075-T6 

76.0 

82.5 

11.4 

4130 

98.5 

117.0 

14.3 


(8) Stress-Strain Curves: Not given. 


Figure 5. - Sample abstract for report containing fatigue data. 


REFERENCE NUMBE R 15 (MCIC 73988) 


Materials: 300 M 

Pendleberry, S. L., Simenz, R. F., and Walker, E. K., "Fracture Toughness and Crack 

Propagation of 300 M Steel", Technical Report No. DS-68-18, Lockheed-California Company, 

August (1968). 

Test Information 

(1) Fracture and Crack Propagation Tests: Tests were conducted on one material in 

three product forms to study the effects of material thickness and strength 
level. 

(2) Type of Test Machine: Lockheed-designed closed-hoop servohydraulic fatigue 

machine (150,000 lb. capacity), Lockheed-designed axial load resonant fatigue 
machine (250,000 lb. capacity), and a universal hydraulic testing machine 
(60,000-400,000 lb. capacity). 

(3) Number of Specimens: 132 specimens were used to obtain both crack propagation 

and fracture data. 

(4) Stress Ratio: R ■ 0.1 or 0.5. 

(5) Test Temperature and Environment: Tests were conducted at room temperature in 

a moist air or salt spray environment, 

(6) Test Frequency: 20 cps to precrack specimens, 

(7) FCP Data: Presented in basic form. 


Specimen Data 


(1) Melting Practice/Heat Treatment: Specimens were normalized at 1700°F/l-l/2 hours, 

air cooled, austinitized at 1600°F/l-l/2 hours, oil quenched, and double-tempered 
at 500°F to 1050®F depending upon strength level desired. 

(2) Ductility: Not given. 

(3) Fabrication Methods: Specimens were machined from 0.125 inch sheet, 0.5 or 

0.75 inch plate or forgings. Precracking was done by axial tens ion- tens ion 
fatigue generated from an EDM slot. 

(4) Surface Finish: Specimens were left as machined. 


(5) Specimen Dimensions: 


Specimen Type 

Th ick- 

Gross 

Net 

Gross 

Net 

Slot 

Slot 


ness, 

Length, 

Length 

, Width, 

Width, 

Length, 

Width, 


inch 

inch 

inch 

inch 

inch 

inch 

inch 

Surface Crack 

.125 

14 

3,5 

4 

2.25 

.08 

.006 

Surface Crack 

.375 

16 

3.5 

5 

2.25 

.08 

.006 

Surface Crack 

.75 

28 

9.0 

12 

4.5 

.08 

.010 

Through Crack 

.125 

15 

--- 

5 

— 

.5 

— 

Through Crack 

.375 

28 

9.0 

12 

5.0 

.5 

.010 

(6) Chemical 

Composition: See 

report for 

analysis of each 

heat of 

material . 



(7) Tensile Properties: See report for original materials' properties, results of 

heat treatment study and properties of control specimens after heat treatment. 

(8) Stress-Strain Curves: Not given. 


Figure 6. — Sample abstract for report containing fatigue-crack 
propagation and fracture data. 
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Each data file always contains at least Cards 1 and 3. Card 2 is an optional 
card which may be necessary to supplement, clarify or expand Card 1 information 
in particular situations. Card 4 is a particular addenda of crack-growth infor- 
mation (i.e., cycle count and crack size) necessary only for fatigue -crack- 
propagation analysis. 

Title or Lead Card (Card 1) Format .— The format of Card 1 is illustrated in 
figure 7. Eleven fields of general descriptive information are presented. Their 
contents are as follows: 


Field 123 4 5 6 7 8910 


/ 











0 rt 









Thickne ss 


£ ° 

~c 






TVS, 

TUS, 

or Diam- 


H 

D 

Reference 


Alloy 

Product 


MN/m^ 

MN/m^ 
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Figure 7, — Format for encoding of title or lead card (card 1). 

(1) The type of data contained in the associated data file is indicated 
in columns 1 through 3 by using an alphanumeric format with three 
coding abbreviations : 

FAT — data from cons tant-amplitude-fatigue tests where the con- 
trolled variable is stress or strain and the dependent variable 
is the total number of cycles to complete failure of the speci- 
men (i.e., the fatigue life). 

FCP — fatigue-crack-propagation data from a constant -amplitude 
stress or strain-cycling test where the size of a fatigue crack 
is monitored as a function of the number of loading cycles. 
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FT — data from a monotonic loading (load or displacement con- 
trolled) test to fracture of a specimen with an initial fatigue 
precrack. 

(2) The dimensional units of the recorded data in the file are identified 
in column 4. A blank denotes International System (SI) of Units; a 
value of 1 indicates English units; and a value of 2 indicates CGS 
units . 

(3) The source reference number is listed in columns 5 through 11 in an 
alphanumeric code of the following format: 

”NNNNNNL", 

where N is a numeric character (0 to 9) and L is an alphabet character 
(A to Z). The numeric code corresponds to the references numbers in 
Appendix A. The suffix letter refers to a specific batch of data from 
the referenced document. This reference number is the same as that 
listed on the succeeding data cards. 

(4) The type of material (e.g., aluminum, steel, and titanium) is described 
in columns 11 through 19 using an alphanumeric format. 

(5) The alloy designation (e.g., SAE 4340, 7075-T651, Ti-6A1-4V) is 
given in columns 20 through 29 in an alphanumeric format. 

(6) The product form (e.g., plate, sheet, bar, forging, and casting) is 
listed in columns 30 through 39 in an alphanumeric format. 

(7) The heat treatment (e.g., Q and T, STA, annealed, normalized) is 
described in columns 40 through 58 in an alphanumeric format. 

(8) The TYS^ MN/m^ (or ksi), is given in columns 59 through 63 in a 
fixed point numeric format. 

(9) The TUS MN/m^ (or ksi), is given in columns 64 through 68 in a 
fixed point numeric format. 

(10) Thickness or diame ter , mm (or in.), of the specimen is listed in 
columns 69 through 74 in a fixed point numeric format. For a 
round specimen where this value represents the diameter, columns 
75 through 80 (item 11 below) will be blank. 

(11) Width , mm (or in.), of the specimen is given in columns 75 
through 80 in a fixed point numeric format. For a round specimen 
these columns are blank and the diameter is given in item 10 above. 
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.Subtitle Card (Card 2) Format . -The subtitle card is an optional card pro- 
vided for particular instances where supplementary information is necessary or 
desirable. This is an open-field card whose format is coded alphanumerically 
and read directly as a subtitle to Card 1 in data listings or tabulations. 

Data Card (Card 3) Format . ~ This card contains the principal test parameters 
and results of each test on which data are collected and compiled. Since the 
types of data may represent either fatigue, fatigue -crack propagation, or frac- 
ture tests, three formats are necessary for this card as detailed in the fol- 
lowing subsections. Where similar test parameters are encountered among the 
types of data, common fields have been assigned to the formats. 

Fatigue (FAT) Data Card Format. “The fatigue data card contains 13 fields 
of information listed in the following formats (see fig. 8): 


Field 1 2 3 4567 8 9 10 U 



Figure 8. — Format for encoding fatigue data card (card 3). 

(1) Specimen identification is listed in columns 1 through 8 using an 
alphanumeric format. 

(2) Maximum stress, MN/m^ (or ksi), or maximum strain is listed in columns 
9 through 14 in a fixed point numeric format. The stress or strain 
option is designated by the Field 5 indicator. 

(3) Stress ratio or strain ratio (ratio of minimum to maximum value) is 
listed in columns 15 through 19 in a fixed point numeric format. The 
stress or strain option is designated by the Field 5 indicator. 

(4) Cyclic frequency , Hz, is listed in columns 20 through 24 in a fixed 
point numeric format. 
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(5) An Indicator is given in column 25 to show whether items 2 or 3 
above are in terms of stress or strain. If this column is blank, 
stress is indicated; and if this column contains an "E", strain 
is indicated. 

(6) The type of notch configuration is listed in columns 26 and 27 
by the following abbreviations: 

CN — center -notched sheet or plate 
EN — edge-notched sheet or plate 
FN — fillet-notched sheet or plate 
CR — circumferentially notched round bar. 

These columns are blank for an unnotched specimen. 

(7) The theoretical stress -concentration factor of the notch geometry 
is given in columns 28 through 32 in a fixed point numeric format. 

These columns are blank for an unnotched specimen. 

(8) The notch root radius , mm (or in.), is given in columns 33 through 37 
in a fixed point numeric format. These columns are blank for an 
unnotched specimen. 

(9) The fatigue life , cycles, is given in columns 38 through 47 in a 
fixed point numeric format. 

(10) An Indicator is given in column 48 to show whether or not the specimen 
was a runout. A "1" in column 48 indicates that the specimen did not 
fail (DNF). 

(11) This is an open field and columns 48 through 69 are left blank. 

(12) Test temperature . °C (or °F), is listed in columns 70 through 73 in 
a fixed point numeric format. 

(13) The source reference number is given in columns 74 through 80 in an 
alphanumeric format of the following type: 

"NNNNNNL", 

where N is a numeric character (0 to 9) and L is an alphabet character 
(A to Z). The numeric code corresponds to the source reference numbers 
in Appendix A. The suffix letter refers to a specific batch of data 
from the referenced document. This source reference number is the same 
as that listed on the corresponding Number 1 Lead Data Card. 

Fatigue-Crack Propagation (FCP) Data Card Format. —The complete recording 
of fatigue-crack-propagation data requires two different card formats. Card 3, 
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described herein, contains the basic test information; Card 4, described later, 
contains the cycle counts and crack size measurements as determined from the 
test. Thus, the data file from a single fatigue -crack-propagation test is 
made up of one Card 3 and one or more Card 4's. 

The layout of Card 3 for the fatigue-crack-propagation test parameters is 
shown in figure 9. 
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Figure 9. ~ Format for encoding fatigue -crack-propagation data card 
(card 3). 

A total of 16 fields are indicated. The field contents are as follows: 

(1) Specimen identification is listed in columns 1 through 8 using an 
alphanumeric format. 

(2) Maximum cyclic stress . MN/m^ (or ksi), or maximum cyclic load, kN 
(or kips) is listed in columns 9 through 14 in a fixed point numeric 
format. The stress or load option is designated by the Field 5 
indicator. 

(3) Stress ratio or load ratio (ratio of minimum to maximum values) is 
listed in columns 15 through 19 in a fixed point numeric format. 

The stress or load option is designated by the Field 5 indicator. 

(^) Cyclic frequency . Hz, is listed in columns 20 through 24 in a fixed 
point numeric format. 

(5) The specimen type is indicated in column 25 as a numeric code and 
supplemented in columns 26 and 27 by an acron 3 mionic code for easier 
identification. Since the specimen type also determines the usual 
convention for selecting either stress or load in the analysis, it is 
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this field that designates the stress or load option for Fields 2 
and 3. The following convention is used in fatigue-crack propa- 
gation: 


Code 

Specimen Type 

Option 

ICT 

Compact Tension 

Load 
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Center Crack 

Stress 
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(6) Specimen thickness , mm (or in.), is listed in columns 28 through 32 
in a fixed point numeric format. 

(7) Specimen width , mm (or in.), is listed in columns 33 through 37 
in a fixed point numeric format. 

(8) Tensile yield strength , MN/m^ (or ksi), representative of that speci- 
men material is listed in columns 38 through 42 in a fixed point 
numeric format. 

(9) Tensile ultimate strength , MN/m^ (or ksi), representative of that 
specimen material is listed in columns 43 through 47 in a fixed 
point numeric format. 

(10) A reference dimension , mm (or in.), is listed in columns 48 through 
52 in a fixed point numeric format. This dimension is utilized when 
experimental measurements are recorded relative to a point other than 
the crack origin prescribed by the analysis. 

(11) This is an open field and columns 53 through 57 are left blank. 

(12) Elastic modulus , 10^ MN/m^ (or 10^ ksi), is listed in columns 58 
through 61 in a fixed point numeric format. 

(13) Specimen compliance , 10“®isr^ (or 10”®lb“^), used specifically for the 
double cantilever specimen is listed in columns 62 through 66 in a 
fixed point numeric format. 

(14) Poisson*s ratio for elastic deformation is listed in columns 67 
through 69 in a fixed point numeric format. 
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(15) Test temperature , °C (or °F), is listed in columns 70 through 73 
in a fixed point numeric format. 

(16) The source reference number is given in columns 74 through 80 in an 
alphanumeric format of the following type: 

"NNNNNNL", 

where N is a numeric character (0 to 9) and L is an alphabet character 
(A to Z). The numeric code corresponds to the source reference 
numbers in Appendix A. The suffix letter refers to a specific batch 
of data from the referenced document. This source reference number is 
the same as that listed on the corresponding Number 1 Lead Data Card. 

Fracture (FT) Data Card Format. ~ Fracture data for a variety of test speci- 
men configurations is accommodated on the card format shown in figure 10. The 
detail presented is dictated to a large degree by the number of important crack 
lengths and stresses which are associated with and reported for thin sheet (plane 
stress) fracture studies. 


field 1 2 3 4 5 6 7 8 9 10 1 1 12 13 14 15 16 


/ 

specimen 

Thick- 
ness, 
mm 
{or in.) 
F5.0 

Width, 
mm 
{or in.) 
F5.0 

Initial 

Crack 

Length, 

Pop- in 
Stress 
MN7m2 
(or ksi) 
or Load 
kN 

{or kips) 
F5.0 

Offset 
Stress 
MN/m2 
(or ksi) 
or Load 
kN 

(or kips) 
F5.0 

Visu- 
ally 
Deter- 
mined 
Critical 
Crack 
Length , 

Photo- 

Record- 

ed 

Critical 

Crack 

Length, 

Maxi- 
mum 
Stress 
MN/m^ 
(or ksi) 
or Load 
kN 

(or kips) 
F5.0 

TVS, 
MN/m^ 
(or ksi) 
F5.0 

TUS, 
MN/m^ 
(or ksi) 
F5.0 

Special 

Dimen- 

sion 

(1), 

e 

1 

41 

Q 

V3 

a 

H 


Test 

Tem- 

pera- 

ture, 

Reference 

Nymber 

A.7 

cation 

A8 

{or in.) 
F5.0 

(or in.) 
F5.0 

(or in.) 
F5.0 

(or in.) 
F5. 0 


< 

Open 

(or F) 

14 

2 4 6 0 

10 12 

14 16 IB 

20 22 

24 26 20 

30 32 

34 36 38 

40 42 

44 46 48 

50 52 

54 56 58 

60 0:2 

04 

66 



^4 -6 

1 1 ! 1 ! j J 

1 M 1 

1 1 1 1 

I ! ! ! 

1 I ! 1 

1 1 M 

1 1 1 1 

. . 1 . 

! ! ! 1 

! ! ! { 

! ! ! ! 

1 1 M 

! 1 ! ! 



~rr 

rT" j T ^ 

! ! ! ! 1 I 


Figure 10. Format for encoding fracture data card (card 3). 

A total of 16 fields of data are contained on the card. Their contents are as 
follows : 

(1) Specimen identification is listed in columns 1 through 8 using an 
alphanumeric format . 

(2) Specimen thickness , mm (or in.), is listed in columns 11 through 13 
using a fixed point numeric format. 

(3) Specimen width , mm (or in.), is listed in columns 14 through 18 
using a fixed point numeric format. 
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(4) Initial crack length , mm (or in.), as measured for the fatigue 
precrack is listed in columns 19 through 23 using a fixed point 
numeric format. 

(5) ' *Pop-in*’ stress , MN/m^ (or ksi), or ' ^pop-in*^ load , kN (or kips), 
is listed in columns 24 through 28 using a fixed point numeric 
format. The stress or load option is designated by the Field 

13 indicator. 

(6) Offset stress . MN/m^ (or ksi), or offset load . kN (or kips), 
is listed in columns 29 through 33 using a fixed point numeric 
format. The stress or load option is designated by the Field 
13 indicator. 

(7) Visually determined critical crack length , mm (or in.), is listed 
in columns 34 through 38 using a fixed point numeric format. 

(8) Photo-recorded critical crack length , mm (or in.), is listed in 
columns 39 through 43 in a fixed point numeric format. 

(9) Maximum stress . MN/m^ (or ksi), or maximum load . kN (or kips), 

is listed in columns 44 through 48 in fixed point numeric format. 
The stress or load option is designated by the Field 13 indicator. 

(10) Tensile yield strength , MN/m^ (or ksi), representative of that 
specimen material is listed in columns 49 through 53 in a fixed 
point numeric format. 

(11) Tensile ultimate strength , MN/m^ (or ksi), representative of that 
specimen material is listed in columns 54 through 58 in a fixed 
point numeric format. 

(12) A special dimension , mm (or in.), characteristic of that specimen 
type is listed in columns 59 through 63 in a fixed point numeric 
format. 

(13) The specimen type is indicated in column 64 as a numeric code and 
supplemented in columns 65 and 66 by an acronymonic code for easier 
identification. Since the specimen type also determines the usual 
convention for selecting either stress or load in the analysis, it 
is this field that designates the option for Fields 5, 6, and 9. 

The following conventions that were used in fatigue-crack propa- 
gation are also used for fracture: 
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Code 

Specimen Type 

Option 

ICT 

Compact Tension 

Load 

2CC 

Center Crack 

Stress 

3SF 

Surface Flaw or Part- 



through Crack 

Stress 

5 DC 

Double Cantilever Beam 

Load 

6NB 

Notch Bend 

Load 


(14) An open field is in columns 67 through 69. 

(15) Test temperature . °C (or °F), is listed in columns 70 through 
73 in an integer format. 

(16) The source reference number is given in columns 74 through 80 in 
an alphanumeric format of the following type; 

"NNNNNNL", 

where N is a numeric character (0 to 9) and L is an alphabet character 
(A to Z). The numeric code corresponds to the source reference numbers 
in Appendix A. The suffix letter will refer to a specific batch of 
data from the referenced document. This source reference number is the 
same as that listed on the corresponding Number 1 Lead Data Card. 

Crack-Growth Card (Card 4) Format .- The crack-growth card is used for record- 
ing the crack-size measurements and cycle counts associated with a given fatigue- 
crack-propagation test or test specimen. Each card contains one set of data 
points. The format of Card 4 for crack-growth measurements is illustrated in 
figure 11. 


Field 12 3 4 



Figure 11. — Format for encoding crack-growth data card (card 4). 
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A total of 4 fields are indicated. Their contents are as follows: 

(1) Specimen identification is listed in columns 1 through 8 using 
an alphanumeric format. 

(2) Number of cycles associated with the first data point on the 
card is listed in columns 11 through 20 in an integer format. 

(3) Crack length , mm (or in.), associated with the first data point 
on the card as measured in the width dimension of the specimen 
is listed in columns 21 through 30 in a fixed point numeric 
format. 

(4) Crack depth , mm (or in.), as measured into the thickness of the 
specimen is listed in columns 31 through 40 in a fixed point 
numeric format. 


Data Retrieval and Sorting 

The data handling system consists of two sets of programs. The first set 
implements the storage of fatigue, fatigue-crack propagation, and fracture data 
on magnetic tape. The second set implements the retrieval of data on the basis 
of certain specified parameters. 

The storage program writes the data in card-image format on seven-track 
magnetic tape at a density of 800 bits per inch. Materials are separated from 
each other by end-of-file cards. There is a different tape for each type of 
data. 

Data retrieval is implemented through a set of programs that sorts the data 
by a number of parameters including stress ratio, stress, frequency, environment, 
and test temperature. These parameters must be specified on a separate control 
card. Specified data then may be transferred from the magnetic tape to any of a 
number of output devices. Information may be obtained in the form of magnetic 
tape, punched cards, or printed output. Either SI or English units may be used. 

Additional analytical subprograms are added to the sorting program to obtain 
graphical output and to perform curve fitting and statistical analysis. Figure 
12 presents a flow chart outlining the data storage and retrieval system. 
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Figure 12. — Data storage, retrieval, and sorting system flow chart. 




STATISTICAL METHODS OF ANALYSIS 


The phenomenological approach to the study of fatigue and fatigue-crack prop- 
agation is usually concerned with formulating a model of material behavior. In 
the present program, this model took the form of a regression equation that was 
fitted to empirical data. Statistical analysis provided a method by which the 
performance of the various empirical models could be compared and evaluated. 

The method which was developed for the analysis of fatigue data is outlined 
in Appendix C. The formulations which were used in the analysis are discussed 
in the following section. 

In the fatigue analysis, a third-order polynomial proved to be useful for 
many of the initial comparative studies. The equation was written in the following 
form: 

log — Aq + "h ^2^eq ^®eq * 

Further investigations revealed that equation (1) could be simplified to a linear 
regression equation involving a single independent variable, 

Y = Ao + A,X , (2) 

where X represented a mapping function linearly related to the dependent variable 
Y. This same simple functional form was also found to be useful in the fatigue- 
crack -propagation analyses.* A least-squares regression procedure was used to 
establish optimum coefficients fur equations (1) and (2). The optimization 
procedure was based on a minimization of the standard error of estimate, 


s 


2(Y-Y)^ 

n-2 


(3) 


Different formulations for the independent variables were compared through 
calculation of the statistical parameter, , This factor, which provided a 
quantitative estimate of goodness of fit, was used to describe the fraction of the 
sum of squares of deviations of the dependent variable from its mean associated 
with the regression. It was defined by the relationship 




1 - SSD/TSS = 1 




(4) 


* The exact definition of this mapping function was omitted here for simplicity, 
it is detailed in the later sections on fatigue and fatigue-crack propagation. 
Briefly, however, for the fatigue analysis X = f(egq) and Y = log Nf, and for 
the fatigue-crack-propagation analysis, X = ^ (Kg££)^and Y = log da/dN. 
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Values of approaching 100 percent were considered most desirable, since 
that tendency indicated a large percentage of the variance of the dependent vari- 
able was attributable to the regression. In the case of fatigue, such values of 
indicated a good correlation between equivalent strain and fatigue life. In 
the case of fatigue-crack propagation, they indicated a well-defined relationship 
between the effective stress-intensity factor and crack-growth rate. The param- 
eter was used extensively in the two sections of this report on fatigue and fatigue- 
crack propagation. 

In parts of the fatigue analysis where the degree of fit for one material was 
compared with that for another material, a modified value of R^, designated as R^^ , 
was used. Such a statistic was necessary when the sample population was fairly 
small in comparison with the number of degrees of freedom. For example, R^^ was 
used to compare a small sample of fatigue data on Ti-6A1-4V alloy with a large 
sample of data on 300M steel, (See Table 2 on page 38.) This term provided a more 
realistic estimate of fit than R^, since it accounted for the number of degrees of 
freedom and the sample population. It provided a sample estimate of the fraction 
of the variance of the dependent variable attributable to regression (ref, 10) and 
was expressed as follows: 


Rm^ - 1 


a-R^)(n-l) 

(n-v-1) 


(5) 


In cases where n was only slightly larger than v, the statistic was appreciably 
smaller than the statistic. However, when n » v, the value of Rj„^ approached 
that of , Thus, only values of were computed for comparison of results from 
large data sample populations. 

After screening the formulations of interest, it was considered desirable to 
establish tolerance limits on the best empirical models. These tolerance limits 
are calculated to define an interval which can be claimed to contain a specified 
proportion of the data population with a specific degree of confidence. Before 
tolerance limits could be calculated, it was necessary to determine whether the 
data satisfied the appropriate statistical conditions. Primarily, the data had to 
be independent and be normally distributed about the regression line and had to 
have zero mean deviations from that line and have a constant variance (ref. 11). 

When the residuals (or deviations from the mean curve) were plotted as a 
function of actual values of the dependent variable, it was possible to determine, 
by inspection, whether the data were independent and had an essentially uniform 
variance throughout their range. Actual values for the dependent variable were 
used since it was then possible to compare different fitting functions without 
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changing the fatigue life values of individual data points on the residual plot. 

If only one fitting function had been studied, it would have been reasonable to use 
predicted values of the dependent variable as is customary in most statistical 
analyses. The additional criterion of log-normality was tested in several cases 
through construction and examination of frequency distribution plots of the 
residuals. Although log-normality of the data was not proved, the frequency 
distribution plots indicated that the data were not skewed appreciably and were 
reasonably approximated by a log-normal distribution. 

After statistical conditions were satisfied, it was possible to calculate the 
estimated variance of a specific value of Y about the regression line. For a first 
order equation [eq. (2)] , the point estimate of variance was defined according to 
the following expression: 


V = s^ 


ri + n 

Ln S(X.X.)2j 
1 


( 6 ) 


To solve equation (6), it was necessary to determine values of X^ for each 
data value based on values of Then X was calculated as a simple average of 

the X^'s. The same process was used in order to calculate s^, based on values of 
After these calculations were completed, the variance was calculated for a 
selected value of X, Knowing the estimates of variance at X, it was then possible 
to determine tolerance limits of level (u) at a desired confidence (v) , according 
to the following formulation (ref, 12) 


n > — 

Y = Y + t Js^+V 

u, V - V 


(7) 


Equation (7) was only valid, however, for data sets which were essentially of 
uniform variance throughout the range of Y* 

In cases where the variance was nonuniform, it was necessary to modify the 
residuals through the use of a weighting function, W(X), so that the transformed 
residuals were approximately uniform. A discussion of this process is included in 
the Fatigue Analysis section of this report. 

In defining all of the above equations, it was assumed that the data under 
analysis were constituents of a single population, presumably from a single source, 
where factors such as between laboratory and between test machine variability were 
of no importance. In a practical situation, however, a large data accumulation for 
a given material is often the result of work at numerous laboratories. In such a 
case, it is inevitable that some portion of the observed data variance is really 
caused by between-system variations. It is desirable to isolate these two factors 
so that the material scatter can be considered apart from the laboratory-introduced 
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scatter. If, for example, it is found that the be tween- laboratory variance of two 
combined data sets is much larger than the within- laboratory variance of either 
subpopulation, it is reasonable to analyze the data separately since the ratio of 
variances indicates that there is a strong possibility the two materials are 
different or the procedure used to test them was not the same. The following para- 
graphs discuss this problem as it pertains to populations of fatigue data. Basic 
aspects of the discussion are also applicable to fatigue-crack-propagation data. 

A method proposed by Mandel and Paule (ref, 13), involving the interlaboratory 
evaluation of a material with an uneven number of specimens from different sources, 
was considered as a means of properly accounting for within- laboratory and between- 
laboratory variance factors. This method was found to be useful in certain cases 
where fatigue data for a given material, although generated at different labor- 
atories, were obtained from tests run at consistent values of stress ratio and 
notch concentration. It was questionable whether the approach had application for 
most of the accumulated data file, however, since the majority of data from differ- 
ent sources were nonuniform in values of stress ratio and notch concentration. To 
use the method for data such as these, it would necessarily have followed that the 
means of consolidation on R and was sufficiently good that individual sets of 
data could not be statistically isolated. This then implied that data from differ- 
ent sources, even though possibly of nonequal stress ratio or notch concentration, 
could have been compared, after consolidation, as identical data. 

Investigations did not provide sufficient evidence to support this conclusion. 
Most sources contained data at only a few values of R. In some cases, a particular 
stress ratio was represented by only a few nonreplicate tests. The same was true 
for much of the data generated with as a variable. This lack of uniform and 
consistent data made it difficult to conclude with confidence that consolidated 
data run at different R and values were completely homogeneous. 

Despite this problem, it was considered appropriate to calculate tolerance 
limits on the combined data sets according to equation (7), since all requirements 
involving randomness, normality, and uniformity of variance appeared to be met 
satisfactorily. Since it was concluded that subpopulations could not be accepted 
or rejected on the basis of an examination of variances, particular data sets were 
included or excluded on the basis of their overall effect on the quality of fit 
(R^) which was obtainable. In some instances, a visual examination of the plotted 
data was sufficient to exclude a particular data subset. 
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FATIGUE ANALYSIS 


Designers of aircraft structural components usually base their fatigue 

analysis on data from stress versus number of cycles to failure (S-N) curves. 

Data for these S-N curves are obtained from constant-amplitude fatigue tests of 

simple notched or unnotched specimens. The stress value in the S-N curve is 

usually either S or S and the S-N relationship is defined for a constant 
mdx ^ 

value of S or R. Curves are generated at several values of S or R to deter- 
m m 

mine the effect of mean stress or stress ratio. To obtain estimates of fatigue 
life for other values of mean stress, interpolations between existing data must 
be made. Average S-N curves are often used to construct modified Goodman dia- 
grams to aid in making these interpolations. A set of S-N curves is normally 
required for both smooth specimens and several sets of notched specimens with 
different notch concentrations. 

Determination of a meaningful set of average S-N curves for a material may 
require 100 or more specimens. If a statistically based S-N curve is required 
for each condition, this number could easily increase to 500 or more specimens. 

Since such large amounts of fatigue data are not available, even for well- 
characterized materials, it is desirable to have an analytical method for 
combining data from different S-N curves to obtain a single curve containing 
sufficient data to allow the development of a statistically based S-N type relation- 
ship for each material. 

The following sections describe the analytical formulations and approximations 
which were used in the development of the final analytical model. The problem of 
consolidation of data generated at different mean stresses is considered first. 

Three different formulations of equivalent strain are reviewed and compared. Next, 
the consolidation of notched data is considered. Various methods of estimating 
local alternating, mean and maximum stress levels are described and critically 
analyzed. The final step relates to the establishment of a functional relationship 
between equivalent strain and fatigue life. The overall results conclude the 
section. 


Equivalent Strain Concept in Unnotched Specimens 

It has been found in work done at BCL that the effect of mean stress on 
fatigue life can be reasonably accounted for through the use of an equivalent 
stress (or strain). Equivalent stress is defined by an equation relating two 
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terms that uniquely define constant-amplitude loading conditions. One term 
represents the cyclic stress amplitude in terms of either AS or S , while the 

di 

other term defines the mean stress, either directly as S or indirectly as S 

ni ■' max 

in conjunction with . If equivalent strain rather than stress is used, AS 
and S^ are replaced by Ae and e^. 

The following section includes a derivation of the equivalent strain equa- 
tions which were evaluated as a part of this program. The determination of 
strain amplitudes through usage of the cyclic stress-strain curve is also de- 
scribed. Factors influencing the final choice of an equivalent strain formula- 
tion are discussed at the close of the section. 


Formulations of Interest . - Two general formulations of equivalent stress 
have been reviewed. The first involves an additive combination of two stress 
parameters , 


S =AS“+BS 
eq m a 


( 8 ) 


Equation (8) reduces to a form suggested by Stulen (ref. 14) when B, a, and P are 
set equal to unity. 


S 

eq 


k S + S 
m a 


(8a) 


When both coefficients, A and B, and the exponent P are assumed equal to one, 
equation (8) sim.plifies to another form, originally proposed by Topper and Sandor 
(ref. 15), 


S 

eq 


S + S 
m a 


(8b) 


Since equations (8a) and (8b) are applicable only in cases where stress levels 
are nominally elastic, it was necessary to consider a more general formulation. 

To account for inelastic stress-strain behavior, equations (8a) and (8b) were 

* 

modified to define an equivalent strain so that equation (8a) was transformed 
to 


e = e + k S /E 
eq a m 


(8c) 


and equation (8b) was rewritten as 

e 

eq 


e 

a 


+ S “/E 

m 


(8d) 


* For unnotched specimens, local equivalent strain was considered to be the same 
as nominal equivalent strain. For notched specimens, the determination of a 
local equivalent strain was of prime interest. 
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The second general formulation of equivalent stress involves a multiplica- 
tive combination of stress factors, 


S 

eq 


D S ^ 
max 


(9) 


When the parameters D - 1, q/ + 6 =1, and m = a, equation (9) describes a form 
proposed by Walker (ref. 4) 


S 

eq 


(AS)"* (S 

max 


(9a) 


For inelastic stress -strain response, equation (9a) was modified to the follow- 
ing form: 


/ \ f n \ 1 

e „ = (2e ) (S /E) 

eq a max 


(9b) 


Since e^ was required to define each equivalent strain, it became necessary 
to calculate the strain amplitude in cases where it was not measured during the 
test. The following section briefly outlines this calculation procedure. 


Strain -Amplitude Determination by Use of the Cyclic Stress-Strain Curve . - 

In tests performed under strain control, values for e and S (or S ) were 

a m max 

known. However, in load control tests, only values of S (or S ) and S were 

m max a 

known and e^ had to be calculated. Use of the cyclically stable stress-strain 

curve provided a good estimate of e from known values of S . A logarithmic 

a a 

trilinear approximation of the cyclic stress-strain curve was defined as follows 

S = E e , 0 ^ S (1) 
a a a ’ 

'a"'> < =a ^ • d<» 

^a"'- ^a<2) < S, 

Appropriate values of the equation parameters for the investigated materials are 
presented in table 1. Experimental cyclic stress-strain data from the present 
study are detailed in Appendix B. 

A number of different parameter values are indicated for Ti-6A1-4V because 
cyclic as well as monotonic properties for the material vary greatly, depending 
on processing and product form. When the titanium data were analyzed, the set 
of cyclic and monotonic values which appeared to most reasonably represent the 
cyclic and monotonic stress -strain behavior of the material were used. 


Selection of a Method . — Initial investigations showed that all three 
equivalent strain formulations [eqs. (8c), (8d), and (9b)] provided good mean 
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CONSTANTS USED TO DEFINE TRI-LINEAR CYCLIC STRESS -STRAIN CURVES 
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Annealed condition, based on unpublished BCL data. 

Annealed condition, coarse microstructure, from Gamble (data source ref. 90). 
Annealed condition, fine microstructure, from Gamble (data source ref. 90). 
STA condition, from Smith, et al. (ref. 18). 




stress data consolidations. Subsequently, the three methods were analyzed in 
detail to determine which method gave the best overall results. 

As stated previously, the major objective in selecting an equivalent strain 
formulation was to consolidate fatigue test data generated at different stress 
ratios so that all data for a particular material might be treated as one set 
and be represented by a single curve. Examination of the Stulen and Topper- 
Sandor equivalent strain relations [eqs. ( 8 c) and ( 8 d)] reveals that data gener- 
ated at nonzero mean stress values are adjusted by a factor related to the mag- 
nitude of the equation parameters k or a, so that the data more closely repre- 
sent zero mean stress data trends. The best value of k or a is determined by 
the relative influence on fatigue life of mean stress as compared to alternating 
stress. A high value of k or a indicates a large mean stress effect. 

An analogous situation exists for the Walker formulation [eq. (9b)]. Mean 
stress is not present directly in the formulation but it can be easily intro- 
duced because " ^a’ this case, lower values of m imply greater 

effect of mean stress, since lowering the m value increases the exponent on S 

max 

making a change in S (and therefore S ) more important relative to e . 

max m a 

To determine the best values of k, a, and m, a third order polynomial equa- 
tion [eq. ( 1 )] was fit to selected data sets for each definition of equivalent 
strain. Which of these three constants was optimized, depended upon which 
definition of equivalent strain was used. The constant k was used for the Stulen 
method [eq. ( 8 c)], a' was used for the Topper-Sandor method [eq. ( 8 d)] , and m v/as 
used for the Walker method [eq. (9b)]. The polynomial was used because it fit the 
results quite well in the region of available data and because it provided a 
convenient tool for comparison of the degree of data collapse obtainable for each 
equivalent strain equation. It was found, however, that the polynomial behaved 
unrealistically outside the range of data. This did not inhibit its use as a 
comparative tool, but did create some doubt as to the polynomial’s usefulness in 
providing a functional relationship between equivalent strain and fatigue life. 
This problem will be discussed further in a later section. 

In using equation (1), equivalent strain was treated as the dependent vari- 
able (i.e., Y = logarithm of fatigue life was treated as the inde- 
pendent variable (i.e., X = log N^). The equation coefficients (Aq , , A^, and 

A 3 ) were determined by least squares regression. Depending upon which of the 
three methods was used, an optimum value of the material constant (k, o', or m) 
was determined by iteratively conducting the regression analysis until a minimum 
value of the standard error of estimate was obtained. Results for five differ- 
ent sets of smooth-specimen data are summarized in table 2. All data points 
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COMPARISON OF EQUIVALENT STRAIN FORMULATIONS FOR UNNOTCHED 2024-T4 ALUMINUM BAR, 
2024-T3 and 7075-T6 ALUMINIJM SHEET, Ti-6A1-4V BAR, AND 300 M STEEL BILLET 
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Data reported by Titanium Metals Corporation of America (data source ref. 70), bar stock was in the annealed condition. 
Data reported by Bateh, et al. (data source ref. 14). 




were equally weighted in the analysis, and results for specimens that did not 
fail (i.e., runouts) were excluded. 

After performing the regression analysis for each material with each formu- 
lation, a review of the results did not provide sufficient evidence for the se- 
lection of one method of defining equivalent strain in preference to the other 
two. In four out of five data sets, the Walker method was as good or better in 
terms of than the Stulen and Topper-Sandor equations. In the one instance, 
where the regression fit was poorer, the R® value was still a very high 96 per- 
cent. 

The Walker method was attractive for one other major reason — the formula- 
tion used for stress -ratio compensation in the consolidation of fatigue data was 
exactly analogous to the equation found useful in the consolidation of fatigue- 
crack-propagation data obtained from tests at different stress ratios. In a real 
structure, where flaw initiation and propagation both may represent a significant 
percentage of the useful service life, it is expedient to treat both phases as 
two interrelated parts of a single damage process rather than as separate phenom- 
ena. Therefore, an equation such as the Walker equivalent -strain equation, which 
compensates for stress -ratio effects in the same manner for both initiation and 
propagation, appeared to be the most useful method of the three investigated. 

Further investigations were then conducted using the Walker formulation to 
determine the importance of specifying an exact value for m. Since the m value 
provided a compensation on stress ratio, it seemed likely that the optimum value 
of m determined by regression for a given data set was related to the stress- 
ratio values for which it was optimized. This was found to be true in a regres- 
sion analysis performed on 2024-T3 sheet data, in which R = -1.0 data were ex- 
cluded. This screening of the data reduced the m value from 0.41 to 0.39. 
Although the difference was only slight, it did indicate that an exact specifi- 
cation of m for a particular material was somewhat unrealistic. To choose a 
reasonable approximate value, however, it was necessary to determine how much 

deviation from optimum was allowable in the specification of m before drastic 
reductions in R^ would occur. Figure 13 illustrates the results of a study per- 
formed on unnotched 2024-T3 sheet data. The percentage reduction in R^ is plot- 
ted as a function of the deviation from the optimum value of m. The R? value 
was reduced less than one percent for all values of m within 0.07 of the optimum 
value. Deviations in m greater than 0.07 from the optimum caused substanital R^ 
reductions, with large deviations in m (> 0.20) causing reductions in R^ of over 
10 percent. 
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Magnitude of Deviation From Optimum, m, Vaiue in 
Walker Equivalent Strain Formulation 


Figure 13. Sensitivity of data consolidation to variation in m value. 
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After reviewing the data for all the investigated materials, it appeared 
reasonable to attempt usage of a single optimum m value for all the data. Table 
3 indicates the resultant decrease in I? for each material when an B? value of 
0.40 was chosen. The greatest R? reduction occurred with a titanium data sample 
in which an approximately 0.60 percent reduction was observed. Since even this 
reduction was comparatively small, an m value of 0.40 was used in all later anal- 
yses. 


TABLE 3 

REDUCTION IN RESULTING FROM SPECIFICATION OF m 


Material 

Optimum 
m Value 

Optimum , 

percent 

TSS 

SSD for 
ra = 0.40 

Reduction in R® , 
for m = 0.40, 
percent 

2024-T3 Sheet^ 

0.414 

91.55 

167.9 

14,24 

0.03 

7075 -T6 Sheet^ 

0,403 

86.29 

28.43 

3.894 

0.00 

300M Billet^ 

0.366 

77.02 

33.02 

7.725 

0.41 

Ti-6A1-4V Bar*^ 

0.426 

86.09 

13.87 

2.012 

0.59 


Data reported by Grover, et al (data source ref. 1). 

Data reported by Bateh, et al (data source ref. 14). 

Q 

Data reported by Titanium Metals Corporation of America (data source ref. 70), 
bar stock was in the annealed condition. ’ 


Local Stress and Strain Approximations 
in Notched Specimens 

Beyond the consolidation of smooth-specimen data through mean-stress com- 
pensation, it was also of interest to combine notched-specimen data in a similar 

fashion by appropriately accounting for notch effects. Since values of e and 

a 

S were used to calculate equivalent strain values for unnotched specimens, it 
max 
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also seemed reasonable to calculate equivalent strain values for notched speci- 
mens in exactly the same way by using adjusted values of e and S , which 

a max 

would be representative of local strain amplitudes and maximum-stress levels. 


Estimation of Local Alternating Stresses and Strains . — Smooth -specimen 
simulations of local stress-strain behavior in notched specimens (ref. 19) indi- 
cate that combined strain hardening and stress relaxation often occur at the 
notch tip during constant -amplitude nominal-stress cycling. To estimate stable 
local values of alternating stress and strain from nominal values, it is neces- 
sary to compensate for this combined hardening and relaxation. Research (ref. 

20) has shown that the effects of strain hardening or softening can be accounted 
for by using a cyclic stress -strain curve in combination with nominal alternating 
strain values modified by an appropriate notch-concentration factor, such as K^, 
Kj, or Kg. 

All three modifying factors were investigated to determine which one gave 
the most reasonable indication of the local strain concentration. Consolidation 
of notched -specimen fatigue data was considered to be a measure of how well local 
strain was estimated. 

as a Strain-Concentration Factor.— The theoretical stress -concentration 
factor was used extensively in initial investigations as an estimate of the 
effective strain magnification at the notch root. In this way, local strain 
amplitude was estimated as follows: 


which is equivalent to the more general form 

e = K^e 
a t a 


(11b) 


when nominal strains are elastic. This method was found undesirable in further 
investigations because estimated strains were unrealistically high and conserva- 
tive in cases where conditions of high nominal stress amplitude and high 
existed. 

The Neuber Rule and Ke.— As an alternative, the strain-concentration factor 
was next investigated. This factor can be calculated by several different 
methods, but the most commonly used method employs a relationship developed by 
Neuber (ref. 21) which states that 

■'t - • <i2) 
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(13) 


The value of K can be written as 
a 


and can be written similarly as 


o 

a S 

a 


K 

e e 

a 


(14) 


Thus, equation (12) can be rewritten as 


O 6 \% 


/ O 6 \ 

h-itf) 


(15) 


a a 


If equation (11) is used to define and the stress -strain function of equation 

(10) is used to define it is possible to rewrite equation (15) so that e is 

given in terms of known values of K , S , and e , and appropriate values of K, 
rk L a a 1 ^ 

^ 2 ,. , and n^. Three different equation forms may result, depending on whether 

the nominal and local strains are elastic or plastic. 

Case l.-If both local and nominal stress and strain are elastic, equation 
(15) reduces to equation (11a). 

Case 2.— If local stress and strain are plastic and nominal stress and 

elastic, the insertion of the elastic modulus and the stress-strain 
function of equation (10) into equation (15) gives 


(or Kg) "g )+l ^ (SaKt.)^E 

Solving equation (16) in terms of local strain yields 


r 


= exp] ln(Sa^K^.^/EKi(or 


Kg))/(ni(or ng)+l)J 


(16) 


(17) 


Case 3. If both local and nominal stresses are plastic, equation (10) must 
be used for both nominal and local stress-strain behavior, so that equation (15), 
in general form, is rewritten as 


Ki (or Kg 


) ca"' = K, (or K,) Oa"' . 


(18) 


* When equation (10) is used to compute local stress-strain response, and e^ 
are replaced by a and e respectively. 

« cL 
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Simplification of equation (18) reveals that local strain amplitude for the fully 
plastic condition is given by 

Sg = exp[ln(Ki(or '^ 2 )k^s ( oj. Kg))/(ni(or ng) + 1)] . (I 9 ) 

In general, the constant terms, , Kg, n^ , and ng, on the left and right sides 
of equation (19) may have different values and must be treated separately. 

A computerized solution of these equations then yields a means of local 
strain (or stress) determination through application of the Neuber Rule. It 
should be pointed out that Case 3 is rarely encountered in most practical appli- 
cations. It is included in the discussion for the sake of completeness. 

To determine the degree of data consolidation possible using the Neuber 
method of local strain determination, local stresses and strains for notched 
fatigue data were calculated from equations (11a), (17), and (19). Comparison of 
calculated equivalent strains for test data at different values and zero mean 
stress revealed that unrealistically high strain amplitude estimates were calcu- 
lated in cases involving high levels of nominal strain and K^. Since no method 
was found to reasonably account for data at these extreme conditions, this method 
was also considered undesirable, at least when used in the manner outlined herein. 

Kj as a Stress and Strain Concentration Factor. — As a third possibility, use 
of Kf was subsequently tested as a means of local strain estimation. This factor 
can also be written in several different forms; values of Kj calculated in this 
investigation were based on a method proposed by Peterson (ref. 22), 

Kf-l 

• ( 20 ) 

This expression was selected for use because it is simple and has been shown 
(ref. 23) to work reasonably well in comparison to a number of other methods of 
calculating K^. Also, it offered a possible solution to the problems observed 
when using or as a strain multiplier, where data fell further above the 
unnotched equivalent strain curve as the value of increased. 

Analysis of various notched-specimen data sets helped support this idea. 

Using a computer procedure to optimize the value of p for a given material, it 
was possible to account for even the highest values of K^.. Results were good 
enough to warrant the use of this method for determination of local cyclic 
strain amplitudes in all further notched-specimen analyses. 
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Table 4 indicates optimum values of p for the four investigated materials. 
Optimum values for the two aluminum alloys were very similar and it was possi- 
ble to approximate these values at p = 0.18 mm (0.007 in.) with a reduction in 
the optimum value of of less than 0.10 percent. The value of 0.18 mm (0.007 
in.) is considerably below that recommended by Peterson (ref. 22). His value 
of 0.63 mm (0.025 in.) gave R® values almost 2 percent lower than the optimized 
value. This difference was explained in part by the fact that the notched data 
were analyzed independently from the unnotched data. If the unnotched data had 
been included, a higher optimum p value would have resulted, because an increase 
in p would have lowered the overall notched curve, bringing it closer to the 
unnotched curve. It was considered desirable, however, to separate notched and 
unnotched data, since higher p values caused layering of the notched data for 
different K^. 

TABLE 4 

OPTIMUM p VALUES FOR TWO ALUMINUM ALLOYS, A 
HIGH-STRENGTH STEEL, AND A TITANIUM ALLOY^ 


Material 

SSD 

TSS 

R®, 

percent 

— 

Optimum p, 
mm 

(in. ) 

2024-T3 Sheet^ 

5.37 

110.10 

95.13 

0.21 

(0.0083) 

7075 -T6 Sheet’^ 

4.73 

146.30 

96.77 

0.17 

(0.0067) 

300M Forging^ ^ 

16.81 

153.40 

89.04 

0.046 

(0.0018) 

Ti-6A1-4V Bar*^ 

5.10 

21.94 

76.76 

i 

0.020 

(0.0008) 


Data were adjusted according to the weighting function W(X), defined in the 
next part of the Fatigue Analysis section, 
b^ 

Data reported by Grover et al (data source refs. 2 through 4). 
c 

Data reported by Bateh et al (data source ref. 14). 

dr. 

Data reported by Titanium Metals Corp. of America (data source ref. 70) 


The optimum p value for the 300M forging data was 0.046 mm (0.0018 in.). 
This compares with a p value of 0.028 mm (0.0011 in.) which was developed from 
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Peterson's empirical formula (ref. 22) based upon an ultimate tensile strength 
of 2000 MN/m® (290 ksi). Although Peterson's value resulted in only a 0.5 per- 
cent reduction in R® , the optimized value of 0.046 mm (0.0018 in.) was used in 
all 300M data consolidation. Since there was substantial scatter in most of 
the titanium data, a reasonable p value was difficult to define. The optimized 
value of 0.02 mm (0.0008 in.) obtained from a sample of Ti-6A1-4V bar data did 
appear to provide reasonable consolidation on for most data sets, so this 
value was used in all subsequent analyses. 

Estimation of Mean and Maximum Stress Levels. ~ Even after appropriately 
determining a local cyclic strain amplitude, it was necessary to develop a 
method for prediction of the effects of stress relaxation on stable local mean- 
stress values. Smooth specimen, strain-controlled tests performed at BCL (see 
Appendix B) indicated that the stable local mean stress under constant-amplitude 
cycling could be approximated by considering a hypothetical mean stress which 
would develop after an initial loading cycle, if no hardening or softening occur- 
red during that cycle. The development of this method and experimental results 
to determine its validity are presented in the following discussion. 

Two simple methods of predicting local mean stresses were evaluated using 
results from strain-controlled tests in which positive mean strains were main- 
tained. As table 5 indicates, especially for the 7075-T6 sheet, the cyclically 
stable mean stresses were low when Ae was large, but they were higher when Ae 
was small. This reduction of mean stress is related to the amount of plastic 
deformation that occurs in each cycle. 

An understanding of this phenomenon can be found through an examination of 
the material stress -strain behavior under these conditions. Upon initial load- 
ing in tension, deformation will follow the monotonic stress-strain curve (Curve 
A in fig. 14) and a is related to e by some function 

a = f^(e) . (21) 

Deformation upon reversal of the loading direction (Curve B in fig. 14) will be 
influenced by the prior loading. If the influence of previous loading is small, 
Curve B can probably be related to the monotonic stress -strain response. How- 
ever, if this influence is large. Curve B would be more closely approximated by 
the stable cyclic stress-strain curve. For intermediate cases, use of a tran- 
sient cyclic stress -strain curve would be more appropriate. 
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TABLE 5 

COMPARISON OF ACTUAL WITH PREDICTED 
VALUES OF STABLE MEAN STRESS 


Specimen 

Strain 

Ratio 

Total 

Strain 

Range 

Ac tua 1 
Mean Stress, 
MN/m® 
(ksi) 

Predicted 

Mean Stress 

By equation (27a) 
MN/m® 

(ksi) 

By equation (27b) 
MN/m^ 

(ksi) 

2024-T3 Aluminum Sheet 

7 

0 

0.0206 

25 

8.3 

>81. 3 




(3.6) 

(1.2) 

(-11.8) 

8 

0 

0.0153 

7.6 

8.3 

>65 




(1.1) 

(1.2) 

(-9.4) 

6 

0 

0.0101 

15 

20 

5.0 




(2.2) 

(2.9) 

(0.73) 

10 

0.5 

0.0100 

36 

29 

15 




(5.2) 

(4.2) 

(2.2) 

7075 -T6 Aluminum Sheet 

8 

0 

0.0204 

49 

41 

73.8 




(7.1) 

(5.9) 

(10.7) 

9 

0 

0.0152 

43 

58 

106 




(6.3) 

(8.4) 

(15.4) 

7 

0 

0.0101 

160 

188 

198 




(23.2) 

(27.3) 

(28.7) 

11 

0.5 

0.0096 

198 

247 

253 





(28.8) 

(35.9) 

(36.7) 


47 








cr 



Figure 14. ” Schematic illustration of stress-strain 
response during first reversal of loading history. 

Morrow (ref. 24) has pointed out that stable stress -strain behavior after a 
reversal in loading can be approximated by multiplying the cyclic stress-strain 
curve by a factor of 2. Using this observation, the translated stress-strain 
values of Curve B, could be approximated by the monotonic or cyclic stress - 
strain function. 



a'/2 

= f 
m 

(eV2) 


(22a) 

or 

a'/2 

= f 

c 

(eV2) 

• 

(22b) 

It follows from equation (21) 

that 






a 

= f 

(e ) 


(23) 


max 

m 

max 


From equation (23), it can be 

shown 

that 





/ 

c 

= 2f 



(24a) 


max 

m 

max 


or 

/ 

a 

max 

= 2f 

c 

(e' /2) 

max 

• 

(24b) 
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Then, it follows that 



a . = a - a' , 

mm max max 

(25) 

and 

a = (a + 0 . )/2 

m max min 

(26) 

By combining equations 

(24), (25), and (26) one finds that 



a = a -f (e^ / 2 ) 

m max m max ^ 

(27a) 

or 

0=0 - f (e' / 2 ) 

m max c max 

(27b) 


Using equation (10) as the cyclic function (f^) and equation (10) again with 
monotonic parameters as listed in table 6 as the monotonic function (f ), two 

HI 

predicted values [eqs. (27a) and (27b)] of mean stress were computed for each of 
the tests with a mean strain. Results of these calculations are compared with 
the actual stable mean stresses in table 5. Examination of the data shows that 
use of equation (27a) gave the most reasonable predictions for both alloys. 
Equation (27b) gave lower predicted values for 2024-T3 aluminum than did equa- 
tion (27a) because this alloy cyclically hardened. This trend was oppostie for 
7075-T6 aluminum because it cyclically softened. It is interesting to note that 
compressive -mean stresses would be obtained with initial loading in compression. 
Also, it is important to realize that this procedure will not apply to variable- 
amplitude loading because each loading cycle is influenced by the prior cyclic 
history. Thus, a more detailed and complete stress -strain analysis as a function 
of loading history would be required for variable-amplitude conditions. 

Establishment of a Relationship Between 
Equivalent Strain and Fatigue Life 

One of the major goals at the outset of this program was to develop the 
capability to estimate, within a desired confidence, the expected fatigue life 
of a particular alloy, given information on maximum stress, stress ratio, and 
(if notched) notch condition. Toward this end, initial work was centered on 
maximum consolidation of notched and unnotched fatigue data for various combina- 
tions of stress concentration and/or mean stress. The Walker equivalent -strain 
formulation, discussed in earlier sections of this report, was found to be use- 
ful in the consolidation process, and good correlations were established between 
and log through the use of a polynomial expression. 


Sgq = Aq + Aj log + A 2 (log Nj)^ + A 3 (log N^)® 


(28) 
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CONSTANTS USED TO DEFINE TRI -LINEAR MONOTONIC STRESS -STRAIN CURVES 
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All data .^re from references as cited in table 




A major difficulty arose, however, when an attempt was made to establish a 
measure of confidence in the calculated fatigue lives. Since log was the in- 

dependent variable in this equation, it was appropriate to establish limits on 

given log N^, rather than the desired result which would have established 
confidence limits on log N^, given 

Fatigue Life as a Dependent Variable . -In order to eliminate the problem 
discussed above, a variety of new formulations were studied employing log 
as a dependent variable. The desirability of a given formulation was based on 
essentially three factors: (1) predictive capability, (2) simplicity, and (3) 

physical significance. 

The first important aspect, predictive capability, was defined solely on 
the basis of the statistic, which was discussed in the earlier section on 
statistical analysis. Simplicity in a formulation was also an important factor 
since the addition of extra terms in an expression often reduces the signifi- 
cance of the coefficients of original terms. Lastly, the physical significance 
of a particular equation was considered important since a physically meaningful 
equation, in contrast to an empirically derived one, was more likely to be use- 
ful in a general application. This, of course, was true only as long as the 
initial insight was correct and was properly applied. 

The following sections outline a variety of attempted formulations involv- 
ing log Nj as a dependent variable. They briefly summarize the relative merit 
of the various equations as applied to notched and unnotched data used in previ- 
ous evaluations. 


Polynomial Data Fitting .— The first method investigated for establishing fa- 
tigue life as a dependent variable simply involved an interchange of variables 
in equation (28), making an independent variable so that 

log = Ao + A^e^^ + AgS^^s + ^ ^29) 

where represented the Walker formulation as expressed in equation (9b) 

Using equation (29), a regression analysis of selected data sets showed that 
R® values were almost equal to those obtained with equation (28). The problem 
of polynomial uncontrollability outside the range of data still existed, how- 
ever. To partially eliminate this problem, it was thought useful to define 
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an intercept value for the polynomial at % cycle of fatigue life in terms of the 
true fracture ductility, in a manner comparable to that suggested by Morrow (ref. 
24). 

This operation eliminated one degree of freedom in the polynomial and re- 
sulted in slightly better curve definition, but it did not sufficiently improve 
the overall usefulness of the polynomial as a functional relationship between 
and log to warrant its implementation. 

Multivariable Stepwise Regression Analysis .— Much of the analytical work 
performed in the early phases of this program was centered on usage of the three 
formulations for equivalent strain, but there was also some interest in and 
some effort devoted towards the development of alternate functional forms 
which could give comparable or superior data consolidations. Multivariable 
stepwise regression was used to test and compare a variety of factors in regard 
to their usefulness as components of a fatigue life prediction equation. 

Two basic equation forms were reviewed, the first of which involved combi- 
nations of e and S , as follows: 
a m 

log = f(e^, S^) . (30) 

Stepwise multiple regression of equation (30) provided an optimum solution of 
the form 

log Nj = Aq + A^e^^S^E + Age^ + Age^S^E . (31) 

The independent variables are listed in order of significance, with the combina- 
tion e /E providing the most significant increase in and the terms e and 
am a 

e^S^/E providing lesser, yet significant, improvements in R^. Including all 
three variables, the accumulated R^ for the equation using unnotched 7075-T6 
data was 82.0 percent. This was a much poorer consolidation than that obtained 
using equation (29). 

The second general equation was defined so that combinations of maximum 
stress and stress ratio or total strain range could be examined, 

log R or Ae) . (32) 

The optimum solution for this combination of variables was found to be an inverse 
relationship, written as follows: 

Nf = Ao + Ai log + Ag(l+R) • 
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The R2 value for this functional form using unnotched 7075-T6 data was 90.9 per- 
cent which was only slightly less than the value obtained for the third order 
polynomial using the Walker equivalent strain term. However, this formulation 
did not appear to be entirely satisfactory either, since it was observed that 
the regression fit for certain stress ratios was much better than that at others. 
Also, the data could not be displayed as well graphically since it was necessary 
to consider two variables and R) in each plot, as compared to a single 

variable (®gq) polynomial equation. So, investigations were continued in 

an attempt to discover a simpler formulation which would accurately model the 
consolidated fatigue data trends. 


Th®_ Inverse Hy perbolic Tangent Function . — A variety of functions were re- 
viewed in the search for a functional relationship which would provide a useful 
empirical model for consolidated fatigue data. Hyperbolic, exponential, and 
power functions were all investigated and found to be unsatisfactory for fitting 
the complete range of available data. However, the inverse hyperbolic tangent 
function provided a reasonable model of fatigue data trends throughout the life 
range of interest, 10 to 10® cycles to failure. 

"^his function, which was chosen and modified specifically to model the 
sigmoidal shape of the fatigue -crack-propagation (da/dN versus curve, was 

also found to provide a useful model of consolidated fatigue data trends. Since 
this relation is derived in detail in the Fatigue -Crack-Propagation Analysis 
Section, it is outlined only briefly here in terms of its application to fatigue 
data. 

To implement its usage, the following functional form was established: 

log = Aq + Aj tanh"^ [$(e^q)] . ( 34 ) 

The scaling function, §(egq), was appropriately defined as 

logCeufie/eeq^ ) 

^(®eq) ~ log(e^/eg) ’ (35) 

where values of and were selected to appropriately bound the complete 
range of data, as illustrated in figure 15. The upper limit, e^, was found to 
be reasonably represented in most cases by the following approximation: 


€ = € 
u eq 



+ 0.0025 

10 


(36) 


53 




Figure 15. — Schematic illustration of regressed inverse hyperbolic 
tangent curve and appropriate functional limits. 

In cases where representative data did not exist at a fatigue life of 10 

QycigSj it was necessary to estimate a value of Attempts to define a value 

for e in terms of true fracture ductility and true fracture strengths resulted 
u 

in values of s which were unacceptably high. Values for the lower limit, e , 
u ^ _ 

were well represented by an equivalent strain value approximately corresponding 

to a fatigue limit. It was determined accordingly as 

e = e 1 - 0.0005 . (37) 

® 'N^ = 10® 

Once again, in cases where representative data did not exist, a reasonable value 

of e was chosen, 
e 

To avoid error in calculation of the inverse hyperbolic tangent function, it 
was necessary to specify values of and which were higher and lower, respec- 
tively, than any calculated value of equivalent strain. 

The Weighting Function. W(X) . - To calculate point estimates of variance for 
establishing tolerance limits on the inverse hyperbolic tangent function, it was 
necessary (in the case of fatigue) to apply a weighting function to the data in 
order to satisfy the statistical requirement of uniform variance. 
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Before this step could be taken, it was important to identify the pattern 
variance# Examination of plotted residuals revealed that the data 
scatter often increased substantially at either extreme of the fatigue life 
range. The increased variance was most evident in long-life fatigue data. It 
was also in this life range that the slope of the fitted curve was greatest (i. 
e., a small change in the independent variable, e , corresponded to a large 
change in the dependent variable, log N^). Since the variance seemed to be re- 
lated in some fashion to the slope of the curve, it appeared desirable to estab- 
lish a weighting function which would modify the residuals (and, therefore, the 
observed data variance) according to the "steepness" of the fitting function. 

Several functions of this type were reviewed. One suggested formulation 
involving a weight factor proportional to the square of the slope of the fitted 
curve, proved to be useful with some modifications. The weighting function was 
expres sable as 


1/W(X) 


d(Y)~ ! ^ 

-d(X)J 


== rd(log Nf)~ j 2 

L d(e^q) J 


where the derivative for the inverse hyperbolic tangent expression [eq. 
found to be 


(38) 

(34)] was 


d(log Nf) ^ -2.0 A, [log(e„/ee)](log 

dC^eq) f [log(e^/e^)]2 - [logCe^Sg/eg^^ )]2} ' ^ ^ 

To make the weighting functioa [eq. (38)] independent of the absolute slope 
of the curve and dependent only on a ratio of slopes at two points along the 
curve, the derivative [eq. (39)] was normalized through division by a minimum 
value of that derivative. For the inverse h3rperbolic tangent function, the mini- 
mum derivative always occurred at the inflection point of the curve which was 
located midway between the function limits. Therefore, at an equivalent strain 
given by 


e 


eq 


e 

u 


+ e 
e 


2 


(40) 


* With log N£ plotted in the customary fashion, along the abscissa, it appears 
that the slope is actually least in the long-life region. It is useful to 
consider, however, that log = f(e^q), rather than the visually implied 
relationship, = f(log Nf). 

** Based on communication with Lars Sjodahl, General Electric Company, Cincinnati 
Ohio, May 8, 1973. ’ 
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the minimum derivative was defined as* 


rddog 

d(e ) 

eq 


-4.0 [log(e^/e^)](log X) 


+ ®e 




[log(®u®e/®eq^]^} 


. (41) 


Then a combination of equations (38), (39), and (41) resulted In a normalized 
weighting function, expressable as 


W(Ceq) = 



«u ®e 

1 / 

2 

/d(log Nf)\ 

®eq 

\ j 


(42) 


The value of W(e^^) was then bounded between 0.0 and 1.0 with values near 1.0 at 
midrange of the fitted curve and values decreasing (according to the square of 
the ratio of slopes) toward 0.0 at the limits of the curve. 

This function was then applied directly to the residuals In the manner shown 
In Appendix H. Characteristically, the function had almost no effect on data 
falling In the midrange portion of the curve. It did, however, substantially 
reduce the relative magnitudes of the residuals near the extremes of the function. 
The overall effect of this weighting operation was an approximately uniform data 
var lance . 


Results of Fatigue Analysis 

Up to this point, the discussion has dealt with the various considerations 
which were Involved In the development of the overall fatigue data consolidation 
and modelling process. The following paragraphs describe the results of those con- 
siderations . 

A fatigue data consolidation and modelling process was developed through 
which a conglomerate set of fatigue test data at various mean stresses and notch 
concentrations could be consolidated Into a single curve and be reasonably de- 
scribed by a simple analytical expression. Also, statistical considerations were 
applied. Incorporating weight factors, so that probability of survival curves 

* In equations (41) and (42) the vertical slash adjacent to the derivative 
designates an evaluation of the derivative at the Indicated point. 


56 



could be constructed about this consolidated data band from which an estimate of 
simple specimen fatigue life for a given material could be obtained from a single 
plot, once the controlling parameters Ae, K^, and r) were specified. An 

outline in Appendix C provides a step-by-step illustration of this procedure. 

This process was successfully applied to 2024 and 7075 aluminum alloys in sev- 
eral different product forms and tempers and to 300M steel in the forged condition. 
It was also used with marginal success on a Ti-6A1-4V alloy, consisting of numerous 
product forms and heat-treatment conditions. 

In these analyses, notched and unnotched specimen data were treated separately 
because there was a sufficient amount of each type of data to consider them on a 
statistical basis. When the two types were combined, the values were decreased 
by amounts up to about 10 percent. Thus, it would be acceptable to combine notched 
and unnotched results when there are not enough data to analyze them separately on 
a statistical basis. A better correlation of notched and unnotched data would have 
been obtained if a more realistic analysis of notch root stress-strain behavior had 
been available. 

Table 7 summarizes the results of the analyses that were made for each material 
using the final model incorporating the hyperbolic tangent function. Weighted 
values are presented for each combined data set. Also, optimum equation coeffi- 
cients are listed, along with the function limits which were employed in each 
consolidation process. The data source references for each material are included 
in the final column. Graphical displays of the consolidated fatigue data are 
presented in Appendix D. The best-fit regression curve is drawn through the data 
and 90 and 99 percent tolerance curves (95 percent confidence) are drawn below 
that line. Comments concerning individual plots are presented in the introductory 
comments of Appendix D. 


FATIGUE -CRACK-PROPAGATION ANALYSIS 

The determination of the safe life of an aerospace structure must be based 
on a detailed knowledge of the entire continuum of damage mechanisms. This be- 
gins with an understanding of the process of fatigue and its role in leading to 
the initiation of macrocracks, and continues as these macrocracks grow to a size 
which may be critical for the complete fracture of a structure or structural 
component. Once a macrocrack has been initiated, crack growth from the initia- 
tion site, due to continuing fatigue damage, must be predicted in a rational 
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RESULTS OF iTOTCHED AND UNNOTCHED FATIGUE DATA CONSOLIDATION USING THE INVERSE HYPERBOLIC 
TANGENT FUNCTION WITH WEIGHTED EQUIVALENT STRAIN DATA 
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Monotonic and cyclic stress-strain calculations were based on data from Ti-6A1-4V Hot Rolled Bar (see tables 1 and 6). 

^Monotonic and cyclic s tress -s train calculations were based on data from Ti-6A1-4V Plato (see tables 1 and 6). 

'Monotonic and cyclic s tress -s train calculations were based on data from Ti-6A1-4V Cylindrical Forging (see tables 1 and 6). 

Monotonic and cyclic s tress -s train calculations were based on data from Ti-6A1-4V Bar (see tables 1 and 6). 




manner. Since an accurate physical model of crack-tip damage accumulation does 
not exist, a fatigue-crack-propagacion model that accurately characterizes the 
mechanical behavior of the material is generally used. Then the model, which 
summarizes or characterizes experimental results, must be inverted to yield pre- 
dictions of structure life under given loading conditions. 

This section of the report describes the formulation of a phenomenological, 
fatigue -crack-propagation model. The initial subsection on mechanical behavior, 
which describes the basic characteristics of the crack-growth process, is fol- 
lowed by a discussion of the problem of modelling this process. This latter 
subsection is broken down into descriptions of formulating the dependent 
variable, independent variable, and analytical model. Methods of evaluating crack- 
growth rate as the dependent variable are discussed in detail. The independent 
variable portion explains various functions used to account for the effects of 
stress ratio. Based on the results of the above work, various analytical models, 
including a nonlinear analytical expression, are examined in Formulation of an 
Analytical Model. The final subsection details the application of the analytical 
model to five sets of data. 

Observed Mechanical Behavior 

Studies have been conducted in numerous laboratories to obtain fatigue- 
crack-propagation data for various materials. Extensive tests have been per- 
formed by various investigators utilizing center -cracked, compact-tension, and 
surface -flaw specimens. Data have been generated on high-strength steel, alumi- 
num, and titanium alloys under both constant- and variable -amplitude loading 
conditions. The present investigation, however, is concerned only with constant- 
amplitude results. Fatigue-crack-propagation data, recorded in the form of 
crack-length measurements and cycle counts (a^, N^) are not directly useful for 
design purposes since a variety of stress levels, stress ratios, initial crack 
conditions, and environmental conditions are encountered. To make use of these 
data, a fatigue -crack-propagation model must account for the effects of these 
parameters on crack growth and, hence, on specimen life. 

In general, the relationship between crack size and number of applied cycles 
can be represented as a crack-growth curve drawn through the raw data points as 
shown in figure 16. The resulting monotonic curve is described in terms of two 
intervals connected by a transition region. 
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Experimental measurement 
Interpolated curve 
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(1) A region (Region 1) of slow growth over a wide range of cycles 
where the slope is relatively small 

(2) A region (Region 2) of exceedingly rapid growth until failure 
or test termination. 


Families of curves for a given material are generated when the maximum stress, 
stress ratio, or environmental conditions are varied. 

Historically, it has been found convenient to model this crack -damage behav- 
ior as a rate process and to formulate a dependent variable based on the slope of 
the growth curve. The instantaneous rate of change of crack length, or an 
approximation to it, 

Y ^ da ^ Aa 

^ “ dN “ AN ’ 


was chosen as the dependent variable for the formulation of a fatigue -crack- 
propagation model. The independent variable for the process was selected to 
account for the more basic mechanical variables of cyclic stress, stress ratio, 
and crack size. 

An appeal to the theory of linear elasticity has suggested that the damage 
severity at the crack tip might be represented by a stress-intensity factor 
which, in general form, may be written as 

K = S/i* g(a,w) , (44) 

where g(a,w) is a geometric scaling function dependent on crack and specimen 
geometry. As a result, the independent variable is usually considered as some 
function of K and stress ratio, or as originally suggested by Paris et al (ref. 
25), some function of 

AK = (l-R)K . (45) 

If the slope of the crack-growth curve is calculated at the various data 
points, and if the s tress -intensity factor is calculated at these same points, 
then the locus of points can be plotted. These variables are generally 

plotted on log-log scales to obtain a crack -growth -rate curve such as that shown 
in figure 17. Examination of this curve suggests several factors of importance 
that will have to be accounted for in the formulation of a crack-growth model. 

In most materials, there is an upper limit to the crack severity and associ- 
ated critical stress-intensity factor which a material can sustain. At this 
critical value, the crack will propagate unstably. For the rate diagram of da/dN 
versus AK or K , the K value is the terminal (or upper) limit on the abscissa 

m3X o 

as illustrated schematically in figure 17. On a Kjnax basis, the rate of crack 
growth becomes very large as K approaches K , such that the growth-rate curve 
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becomes asymptotic to this limit. On a ^ basis, this limit is (l-R)K . At the 
other extreme, a minimum crack-growth rate of zero is anticipated at a zero value 
of AK or K value. However, this assumption appears to be conservative because 
of evidence that there actually exists a threshold below which there is no 
fatigue-crack propagation. (See fig. 17.) 

The doubly logarithmic plot of da/dN versus K reveals a curve having sig- 
moidal shape. As an approximation, the curve might be represented by three 
linear segments. The first segment, beginning with crack initiation at Kq , is 
steeply sloped and indicates rapid rate of change of crack-growth rate. The 
second segment represents a longer interval of slower rate of change of crack- 
growth rate. The third segment also has a high slope and represents rapid, 
terminal crack growth near K^. Most of the available test data lies in the 
second interval. 

Within the general curve shape, described above, systematic variations in 
the data point locations are observed. When data from tests conducted at sev- 
eral different stress ratios are present, the plot of crack-growth rate versus 
stress -intensity factor will be layered into distinct bands about the locus of 
points having zero stress ratio. (Refer to fig. 3 on page 7.) Layering of data 
points also occurs as a result of variation in such parameters as test frequency, 
environment, and specimen grain direction. It is desirable to predict the char- 
acteristic effects of each parameter; thus, many researchers have formulated 
mathematical models accounting for these parametric effects. Assuming the vari- 
ables K, R, and da/dN, the general form for the fatigue -crack-propagation model 
was established as 

S = f(K, R) . (46) 

The following discussion describes efforts to obtain a useful functional form 
for f(K, R). 


Structure of the Modelling Problem 

The basic concepts discussed in the previous section suggest that the model- 
ling procedure can be thought of as consisting of three distinct steps. 

(1) Formulation of a dependent variable . — How can the crack- 
growth rates be best calculated from the discrete (a^, N.) 
data points? 


63 


(2) Formulation of an independent variable . — What combination 
of R and K can best be used to formulate an independent 
variable that will consolidate the crack -growth-rate data? 

(3) Formulation of an analytical model . ~ What functional form 
containing the dependent and independent variables should be 
chosen to best approximate the sigmoidal character of the 
crack-growth-rate curve? 

An approach to the solution of these three modelling problems is described in 
the following subsections. 

Particular effort was devoted to obtaining an expression in which compensa- 
tions for the effects of stress ratio were uncoupled from the factors influenc- 
ing general curve shape. Such a feature permits a greater flexibility in the 
analysis of fatigue -crack-propagation data. Although several parameters affect 
the distribution of data points, stress ratio is the most significant of these. 
Accounting for other important parameters, such as frequency, was beyond the 
scope of the present work. 

Formulation of a Dependent Variable . — It is necessary to obtain values for 
the dependent variable, the crack-growth rate, from the (a^, N^) data. Two basic 
methods of deriving the crack -growth rate have been used in the past; curve fit- 
ting and incremental-slope approximation. Curve fitting implies that an analyti- 
cal expression is fitted to all of the crack-growth data by least-squares regres- 
sion. Incremental- slope calculation implies the use of a divided differences 
scheme to find the slope at any given point along the crack -growth curve. 

From the analyses conducted, it is apparent that the determination of such 
a derivative by means of some analytical expressions is far less desirable than 
the use of a local or segmental fit to the data. Since this observation has 
been made in all of the data sets analyzed, a formalized illustration of the 
inadequacies of the fitting of a single analytical expression is presented. 

Of the general analytical expressions which are available for curve fitting, 
the most popular choice of functions with respect to numerical considerations are 
polynomials. To explore the application of polynomials in fitting the crack- 
growth curve, one must consider the characteristics of the crack-growth curve. 
Typically, two regions of the curve from crack initiation to specimen failure may 
be described as done earlier and shown in figure 16. These curve segments are 
connected by a transition region having a considerably smaller radius of curva- 
ture. It is observed that over the entire range of cyclic values, the curve is 
monotonically increasing. 
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Some general observations can be made about polynomials that are relevant 
to this situation. Consider a polynomial of degree q, 

a(N) = Aq + AjN + AgN® ... + A^N^ . (47) 

The first derivative of equation (47) is expressed by the following relation: 

da/dN = A^ + 2AgN + ... + (q-l)A^N^‘^ . (48) 

Equation (48) possesses q roots implying a finite number, q - 1, of extrema over 
the range of the function. Since the polynomial is not a strictly monotone func- 
tion over its range, it is necessary to utilize the function on regions where it 
does display monotone behavior. The existence of extrema in a candidate curve- 
fitting function presents very real difficulties. Sections of the curve having 
negative slope due to extrema would represent the physically impossible situation 

of negative crack growth. Such a result is unacceptable. Furthermore, since it 
is generally desirable to obtain the logarithms of the crack-growth rates, log 
(da/dN), also will be undefined at points having negative slope. 

Figure 18 represents attempts to fit second through seventh degree polyno- 
mials to a typical crack-growth curve. These particular data were obtained from 
a 9. 6 -inch -wide, centrally cracked panel of 0. 29-inch- thick, mill-annealed 
Ti-6A1-4V plate tested at a maximum cyclic stress of 5 ksi and a stress ratio of 
0.1. No terms greater than degree seven were added because of computational 
difficulties encountered in dealing with the coefficients. Successive addition 
of higher order terms improved the fit of the polynomial to the data as indicated 
by the sum of squares of deviation listed in table 8. Although the higher order 
terms improved the fit, they introduced extrema with their resulting negative slopes. 
Because of these extrema, the fitted functions obtained are unsatisfactory models 
of crack growth. 


TABLE 8. 

COMPARISON OF CURVE FITTING RESULTS FOR 
POLYNOMIALS OF DEGREE q = 2 TO q = 7 


Degree of Polynomial 

SSD 

2 

7.4038 

3 

3.8250 

4 

2.1020 

5 

1.2387 

6 

0.7809 

7 

0.5254 


Examination of figure 18 reveals a very close fit of data in the region of 


high crack-growth rate for q = 6 and q = 7. In this situation, where only 
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isolated segments of the crack-growth curve are to be dealt with, the polynomials 
gave good results. However, this type of function seems to be unsuited to appli- 
cation to the whole cyclic range. Thus, either the selection of another class of 


candidate functions for curve fitting or an appeal to point-to-point methods of 


slope evaluation was necessary. 

Several other curve-fitting functions were examined, 
exponential series, 

^ iN 

a(N) = S A.e 

i=o ^ 


Among these was an 


(49) 


Cases of q - 1 and q = 3 were tried, but no improvement was found. 

The point-to-point method of slope evaluation involved use of the various 

divided difference schemes. The first divided difference was merely the slope 
between two adjacent points; thus 


DD 


l(i) 


Aa 

AN 


i+1 


^i+l 



(50) 


Since the crack growth achieved between two data points was usually small, 
observational errors, measurement errors, and subtle material variations influ- 
enced the rate evaluation. With this technique, each rate or slope determination 
was defined entirely by the local conditions. An averaging of these variations 
was achieved by using higher order divided-difference schemes. In the next 
level of refinement, the three-point divided-difference technique, a selection 
of successive subsets of three data points was used to specify the derivative 
at the central point. By using Newton's interpolation formula to define a 
second-degree polynomial through the three data points, we may express the 
derivatives at the intermediate point i as 


ft%. Ni+il - NJ 

da/dN|j - + (Ni-Nj.i) 


H-1 


(51) 


where 




and 


a.+l - a. 

"m' ■ N^Vi - Ni 


are the first divided differences. From a physical perspective, this can also 
be viewed as a slope-averaging technique since the first divided differences are 
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merely the slopes between data points. The mechanics of selecting subsets of 
data are shown in figure 19a. Figure 19b presents an array of divided differ- 
ences in which the progression to higher order approximation can be seen. 

A comparison of regression results for 2-, 3-, and 5-point subsets showed 
the superiority of the five-point divided-difference method in evaluating 
fatigue -crack-growth rates. Some researchers have utilized a seven-point di- 
vided-difference technique. Even though this method may result in a slightly 
better evaluation of crack-growth rates, it is questionable as to whether the 
magnitude of improvement would justify the added computational complexity. Use 
of a divided-difference technique implied that a certain number of data points 
had to both proceed and follow the data point at which the slope was being eval- 
uated. Consequently, q - 1 data points had to be discarded when a qth order 
divided difference was used. In data sets where a small number of readings was 
taken, this feature often caused rejection of the entire set. 

Most of the analyses performed in this study, involved use of the five-point 
divided -difference method. 


Formulation o f the Independent Variable . — It was previously suggested that 
the independent variable be some function of K and R. As a general form for the 
independent variable, assume 

• (52) 

where U(R) is a functional relation that accounts for the effect of stress 
ratio. 

A number of different forms for U(R) have been proposed. The simplest of 
these is U(R) = 1.0. In this way, it is asserted that no stress ratio effects 
are present; then, 

^eff ^max * (53) 


This relation is appropriate if no variation in stress ratio is contained in 
the data, or if the material is insensitive to changes in stress ratio. 

The stress -intensity range also may be used as an independent variable. 
Letting U(R) = (1-R), the expression 


K = (l-R)K = AK 
ef f max 


( 54 ) 


results. This relation has been widely used in the past. 
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Walker (ref. 4) proposed that the Independent variable should represent a 
combination of maximum stress intensity and stress -intensity range. Letting 
U(R) = (1-R)“, has the form 


K = (l-R)"* K 
eff max 


(55) 


Mukherjee and Burns (ref. 26), and Roberts and Erdogan (ref. 27) have proposed 
similar relations in their respective studies. 

More recently, Elber (ref. 28) proposed a fatigue -crack-propagation model 
that is based on crack-closure concepts. Elber observed that a crack in a 
center-cracked panel tended to close before the tensile load was removed. As a 
result, he defined a crack-closure stress below which the crack would be totally 
closed. A general form for the crack-propagation independent variable, based on 
these considerations, may be obtained. Elber proposed that U(R) = (1-R)(1+MR), 
so that 

K = K (1-R)(14-MR) , (56) 

eff max \ / 


where M is determined by optimization or by an experimental procedure. 

The four candidates for U(R) may be compared graphically. Since U(R) rep- 
resents a shifting factor accounting for the effect of stress ratio, it is 
reasonable to plot U(R) versus R for the four candidate functions (fig. 20). 

Nominal coefficient values have been chosen in both the Walker and Elber 
relations to represent application to 7075-T6 aluminum alloy data. When U(R) = 
1.0, no shifting for stress ratio occurs. If U(R) = (1-R), then a linearly 
varying shifting factor from U(R) = 2.0 to U(R) = 0 is generated. Setting 
U(R) = (1-R)”^ produces much greater variation in U(R) for positive stress ratio 
than for negative stress ratios. A similar observation is made when U(R) = 
(1+MR)(1-R). 

The selection of a form for the independent variable should be based on 
physical insights as well as on statistical performance. Although physical 
arguments are not completely formulated at this time, some general considera- 
tions are possible. Since it has been observed that most materials exhibit 
stress -ratio dependent behavior, it is reasonable to assume that the choice of 
U(R) =1.0 would seldom be satisfactory. It is also not reasonable to assume 
that U(R) = (1-R), i.e., that the behavior is governed only by the stress- 
intensity range. The Walker formulation, which is a combination of these two 
effects, is a more physically justifiable selection. Taking a rather different 
approach, Elber based his expression directly on the observed physical behavior 
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Figure 20. Comparison of four candidate effective stress-intensity 
functions . 
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of crack closure. Table 9 contains a statistical comparison of the formula- 
tions. A three-point divided-difference scheme was used to evaluate the crack- 
growth rates. The four independent variable forms were used in a linear fa- 
tigue -crack -growth model. It was found on the basis of R^, that the Walker 
expression, followed by the Elber expression, provided the best consolidation 
of data. The Walker formulation was chosen for use as an independent variable 
in the fatigue-crack-propagation model. 


TABLE 9 

COMPARISON OF INDEPENDENT VARIABLE FORMULATIONS FOR REGRESSION 
ANALYSIS OF DATA ON 2024 -T3 ALUMINUM ALLOY 


Formulation 

SSD 

R2 

^eff ” ^max 

257.40 

0.793 

= (l-R) K 

max 

205.99 

0.769 

B 

1 

T— 1 

II 

104.29 

0.915 

= (l-R)(l-fMR) K 

max 

186.80 

0.890 


Since the Walker formulation of the independent variable was chosen for use 
in the following analyses, it is of interest to examine the nature of the coeffi- 
cient m. Investigation of crack-growth-rate curves indicate the formation of 
bands of data with respect to stress ratio as indicated in figure 3 on page 7. 
From a graphical point of view, the coefficient m, which affects the coupling 

between K and AK, caused a shift of the data bands, i.e., a collapse of data 
max 

towards the mean curve. In the case of a set of data having both positive and 
negative stress ratios, the points collapsed toward the R = 0 data since these 
data are not affected by coupling through m. When the subset consisted of two 
positive stress ratios, the coefficient m was selected to produce the best col- 
lapse of the two stress ratios towards a central line between them. 

The value of the parameter m was obtained for various sets of data by mini- 
mizing the SSD value. A series of investigations was undertaken to determine 
the variations of m with respect to stress ratio within a set of data for a 
particular material. The data sets used previously (7075-T6 and 2024-T3 alumi- 
num) were partitioned in various ways for analysis. The results are presented 
in table 10. 

These correlations indicated that the coefficient m was highly dependent 
on the stress ratio distribution and the number of data points. In other words. 
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different m values were obtained when different subsets of data of a given mate- 
rial were regressed. Generally, ^ ^ 

• The formulations da/dN = C(K and da/dN = C(AK)^ were equally 

max 

satisfactory when R = constant. 

• When the subset of data consisted of specimens for which R > 0, m 
tended to be greater than 0.50. 

• When the subset of data consisted of specimens for which R < 0, 
m tended to be less than 0.50. 

A dependence of m on the material properties probably also exists. To un- 
cover this relation, it would be necessary to compare test results for different 
materials. The comparison sets would have to consist of an identical number of 
data points run at the same stress ratio. Unfortunately, data meeting these 
requirements were not available. 


TABLE 10 

VARIATION OF COEFFICIENT m WITH RESPECT TO DISTRIBUTION OF STRESS RATIOS 


Material 

Stress Ratio 

m 

SSD 


7075 -T6^ 

-1.0 to 0.80 

0.37 

38.90 

0.908 

7075 -T6 

0.0 to 0.80 

0.53 

23.70 

0.914 

7075-T6 

-1.0 to 0.0 

0.04 

4. 10 

0.972 

7075 -T6 

0.0 to 0.33 

0.50 

7.30 

0.922 

7075 -T6 

0.50 to 0.80 

0.70 

12.20 

0.912 

7075 -T6 

-1.0 to -0.80 

0.32 

1.77 

0.974 

2024 -T3^ 

-1.0 to 0.70 

0.42 

97.00 

0.920 

2024-T3 

0.0 to 0.70 

0.44 

104.00 

0.915 


^Data from Hudson (data source ref. 48), and Dubensky (data source ref. 32). 


Formulation of an Analytical Model . — Numerous models of the type illus- 
trated by equation (46) have been formulated by researchers during the last 
decade. Excellent reviews of the literature have been presented in papers by 
Erdogan (ref. 29), Hoskin (ref. 30), and Coffin (ref. 31). Several fatigue- 
crack -propagation laws that have been widely used are described below. All of 
these are empirical equations relying upon regression analysis to calculate 
empirical coefficients. These relations can be quite logically divided into 
classes of linear and nonlinear functions. 
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Linear models, of necessity, neglect initial and terminal behavior. The 
general form for the linear model is 

3S ■ • (57) 

Best known of these models is the linear law of Paris (ref* 32), 

|| = C[(l-R)K]'^ = C(AK)"' . (58) 

This equation is commonly fitted to the data in log-log form to yield the Paris 
regression coefficients C and n'. The Paris model is a linear approximation to 
the rate curve that incorporates a term to account for the effect of stress 
ratio. Although the law generally fits only the central segment of the data 
accurately, it has been used extensively in the literature. 

Other linear models are possible, and several have been proposed. These 
relations, which must be considered elaborations on the Paris model, are due to 
Elber (ref. 28), Walker (ref. 4), and Roberts and Erdogan (ref. 27). The pri- 
mary differences in these expressions lie in the choice of the independent 
variable as discussed in the previous section. 

Modifications of the linear Paris model have been made to create a non- 
linearity at the terminal end of the curve. To approximate the sigmoidal char- 
acter of the rate curve, and to better account for the effects of stress ratio, 
Forman (ref. 33) proposed the relation, 

da ^ C[(l-R)Kl*^ 

dN (1-R)(K^-K) . (59) 

Forman’s equation contains a singular term in the denominator to model the 
terminal region of crack growth. As K approaches the critical stress intensity, 
the denominator goes to zero. Manipulation of the Forman equation leads to 

-R)" \ 

K^-K / . (60) 

The term, (1-R)^ is clearly similar to the Walker formulation for the inde- 
pendent variable and as such helps to account for the effect of stress ratio. 
Forman’s equation has no provision for modelling the interval of crack initia- 
tion and, hence, generates only half of the sigmoidal curve. Variations on the 
form of the singularity are possible. 

A computer program was written to evaluate various fatigue-crack-propaga- 
tion laws. This program computed K values and calculated crack-growth rates by 
three-point divided differences. The models were fitted to these results by 
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linear regression. Sets of data for 7075-T6 and 2024-T3 aluminum alloys and 
T1-6A1-4V titanium alloy were used for comparison purposes. Results of these 
regression analyses, in terms of R® values, are compared in table 11. 

Variations on the form of the singularity in the Forman equation proved to 
be ineffective. The linear model with the Walker formulation for the independ- 
ent variable and the Forman model showed the most promising results. 


The other approach to the modelling of the crack-propagation process is to 
assume a nonlinear function. Recently, Collipriest (ref. 34) suggested a fa- 
tigue -crack-propagation law to model the entire rate curve. This nonlinear 
equation is based on the inverse hyperbolic -tangent function. The model may be 
written as 


da _ 


In K -In 


• arc tanh 


In AK - In Kc(l-R) + In AK 


In K^(l-R) - In AKq 


+ In -fc 


exp 


AKp . ^ 


(61) 


In Collipriest *s equation, the independent variable takes the form of K = AK. 

eff 

This nonlinear approach was investigated further because it seemed to provide a 
realistic method for analysis of fatigue-crack-propagation data. Rather than 
utilizing Collipriest * s equation, it was decided to derive a fatigue -crack- 
propagation model that would allow the implementation of the most effective of 
the independent variable formulations described earlier. The goal of this deri- 
vation was also to obtain a more compact analytical form for the fatigue-crack- 
propagation model. 

The model was based on the inverse hyperbolic tangent suggested by 
Collipriest. The shape of the inverse hyperbolic- tangent function is shown in 
figure 21. The functional form assumed was 


log = Cl + C 2 tanh 1 [$(Kgfj) | . (62) 

The coefficients, Ci and Cg , were to be determined by least squares regression. 

Examination of the tanh~^ curve suggested the proper form for f(K ). The 

eff 

function, was chosen to scale values of the effective stress -intensity 

factor into values of the argument, thus positioning the tanh"’- curve relative 
to the rate curve. Figure 21 shows that the tanh"’- function goes to infinity 
at the values of § = -1 and I = +1. The initial and final conditions of the 
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COMPARISON OF FATIGUE-CRACK PROPAGATION MODELS 
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physical problem implied that the rate of change of da/dN should go to infinity 
at both the terminal and threshold values of K. Clearly, the regions of rapid 
rate of change of da/dN should correspond to arguments in the neighborhood of 
± 1. To establish this correspondence, a function was assumed to scale the K 

err 



Figure 21. - Inverse hyperbolic tangent function. 

The physical initial and final conditions were assigned to the points (log 
Kq , -1) and (log K^, +1) on the $ - log K plane as illustrated in figure 22. 

Assuming a linear scaling function, 

i = M'log K + , (63) 

the slope and the intercept were determined by applying the conditions 

$ = 1 when K = log 

$ = -1 when K = log Kq . (64) 

These conditions yielded a system of simultaneous, linear, algebraic equations, 

1 = M'(log K^) + 

-1 = M'dog Ko) + . (65) 
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It should be noted that when data containing only one stress ratio are 
regressed, the formulation provided equally satisfactory results. In this 

case, the model was 



c, + 


tanh”^ 


log [KcKp/Kmax^] ' 
log [KqK^] 


(69) 


Analyses of Data: Application of the Inverse Hyperbolic Tangent Model 


A fatigue -crack -propagation model that successfully accounted for the ef- 
fects of stress ratio made possible the combination of sets of data from differ- 
ent sources. This was particularly desirable since the ultimate goal of the 
modelling effort was the characterization of the crack-propagation behavior of 
specific materials. Accordingly, it was necessary to obtain data over as wide 
a range of stress-intensity factors as was available. This collection effort 
included not only data from different specimen types but also from different 
heats of material. Data from various sources were combined on the basis of 
visual inspection although statistical techniques for combining data sets were 
available. Statistical methods (standard deviation and tolerance limits) were 
applied to regression equations for the combined data sets to complete the mate- 
rial characterization. 

A computer program was developed to apply equation (68) to large sets of 
fatigue -crack-propagation data. Starting with encoded (a., N.) data, this pro- 

r ^da t ■■ 

gram fitted the inverse hyperbolic-tangent model to ^ K values in the 

following steps: 

(1) Crack -propagation rates were evaluated by a five-point 
divided -difference scheme. 

(2) Maximum stress-intensity-factor values were calculated 
using the appropriate formula for the given specimen type. 

(3) Values of the argument were calculated from the 

K results. 

da I 

(4) Equation (67) was fitted to the (^(^^gff)> ^ values by 
least squares regression. The coefficient m was optimized 
by minimizing the SSD value through iterative regression. 

(5) Statistical parameters, including SSD, R^, and S, were 
generated. Tolerance limits were computed from equation 
(7). 

(6) The data, regression mean curve, and 90 and 99 percent 
tolerance-limit curves were plotted. 
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Analysis of the data necessitated that a selection be made for the values 
of and K^. An excellent summary of threshold values is presented in the 
paper by Donahue et al (ref. 35). Data are included in this source on a large 
number of materials. Values of can be found in such publications as the 
Damage Tolerant Design Handbook (ref. 36). These two sources yielded average 

values for the and limits on the crack-growth rate curve. The data sets 

analyzed also contained upper and lower bounds on K 

max 

Nominal values for were selected from the paper by Donahue et al (ref. 
35). Nominal values for K were established by inspection of the K values 
for the data sets. Data on five materials were analyzed: 7075 -T6, 7075-T7351, 

and 2024-T3 aluminum; 300M steel; and Ti-6A1-4V alloy. The composition of these 
five data sets are listed below. 

7075 -T6 Aluminum Alloy. — Data on center-cracked bare and clad specimens 
were compiled from reports authored by Hudson (data source ref. 48), Hudson and 
Hardrath (data source ref. 92), McEvily and Illg (data source ref. 93), Broek 
(data source ref. 118), and Dubensky (data source ref. 32). 

7075-T7351 Aluminum Alloy. — Data for center-cracked bare specimens includ- 
ing wide panels were obtained from unpublished BCL work and from Feddersen (data 
source ref. 41). 

2024-T3 Aluminum Alloy. — A large amount of data on center-cracked bare and 
clad specimens was taken from reports published by Broek (data source refs. 118 
and 119), Hudson and Hardrath (data source ref. 92), McEvily and Illg (data 
source ref. 93), Schijve et al (data source refs. 68, 120, and 121), Dubensky 
(data source ref. 32), and Carter (data source ref. 128). 

Ti-6A1-4V Alloy. — Data on both center-cracked and compact-tension speci- 
mens were extracted from reports by Feddersen (data source ref. 125) and Bucci 
et al (data source ref. 115). 

300M Steel Alloy. - Data on center-cracked specimens tested in humid air and 
saltwater spray environments, which covered a limited range of stress-intensity 
factors, comes from a report by Pendleberry et al (data source ref. 15). 

Detailed results of the regression analyais performed on the five materials 
are presented in Appendix F and summarized in table 12. Number of data points, 
regression and optimization coefficients, and values, and statistical 
parameters are presented. Appendix figures FI through FlO show the consolidated 
fatigue -crack-propagation data, the fitted curve and tolerance limits, and the 
plotted residuals. 
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Good characterizations of the data were obtained in most cases. Particu- 
larly satisfactory results were achieved for the titanium alloy (fig. F9). The 
sigmoidal character of this crack-growth -rate curve is clearly displayed. 

Rather poor results were obtained for the 300M steel (fig. F7). These data 
included only a limited range of stress -intensity factors and contained a large 
amount of scatter. This scatter probably represents the inherent behavior of 
the material because similar observations were made earlier for fatigue data on 
300M steel. 

A final comparison between three methods of fatigue-crack-propagation ana- 
lysis was made. The five data sets were regressed in three different ways; with 
the inverse hyperbolic -tangent model, with the Paris model [eq. (58)], and with 
the Forman model [eq. (59)]. The results of the comparison are presented in 
table 13. From this table it is observed that the inverse hyperbolic-tangent 
model provided significant improvement in representation of the data, compared 
with the other two methods. 


TABLE 13 

COMPARISON OF FATIGUE -CRACK-PROPAGATION 
MODELS WITH COMPLETE DATA SETS 


Values of R® 


Material 

, da 

C + N log K(l-R) 


+ Cg tanh ^ $ 

7075 -T6 


0.669 

0.875 

0.912 

2024-T3 


0.829 

0.877 

0.923 

Ti-6A1-4V 


0.939 

0.970 

0.982 

300M 


0.415 

0.585 

0.661 

7075-T7351 


0.880 

0.926 

0.952 


FRACTURE ANALYSIS 

The accumulation of damage in a structural material terminates at fracture 
instability. In a chronological sense, this event concludes a chain of crack- 
damage processes such as have been portrayed in the previous sections. Fracture 
toughness and residual strength provide the quantitative characterizations of 
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fracture instability which are necessary to determine both the load carrying ca- 
pacity of the material at a given stage of crack damage and the life remaining 
for subcritical damage processes. 

Fracture Toughness and Residual Strength 

Although fracture toughness and residual strength are related descriptors 
of fracture instability, they do imply different subtleties in the fracture event 
as well as different perspectives on the occurrence of crack extension prior to 
fracture. Generally, fracture toughness refers to a distinct material character- 
istic associated with abrupt fracture instability, under a rising load, after 
only minimal amounts of crack extension. In contrast, residual strength refers 
to fracture behavior which is accompanied by much larger amounts of crack exten- 
sion prior to the critical instability. The former term is usually associated 
with relatively brittle fracture under quasi-plane-strain conditions of stress 
state in the material, while the latter is associated with quasi -plane -stress or 
trans itional-s tress -s tate behavior . 

In name, the term residual strength infers that useful strength remains in 
the structural material even after some stable extension of the crack. As will 
be seen later in the discussion, residual strength is also quantified in dimen- 
sions of toughness and is frequently identified as "apparent" fracture toughness. 

Factors Influencing Fracture Behavior 

There are a large number of material, metallurgical, and mechanical vari- 
ables which influence fracture behavior. These include alloy composition, pro- 
cess details (i.e., mechanical reduction and/or heat treatment) associated with 
a product form, the stress-state effects related to product size, and tempera- 
ture. 

While such an itemization of primary factors may suggest that a characteri- 
zation of fracture behavior can be achieved through a simple categorization of 
these details, such is not the case generally. These factors are highly inter- 
dependent, and a discrete segregation of effects is frequently impossible or, 
at least, not economically feasible. 

For example, for a given alloy composition and product form, a specific 

section size (and, hence, stress-state characteristic) may be associated with a 
particular degree of mechanical reduction, such that another section size may 
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have a distinctly different level of mechanical reduction. In other words, two 
different size product forms of common alloy composition may, in reality, be two 
different materials from the perspective of fracture behavior due to differing 
degrees of contained mechanical work. Similarly, in quench-rate sensitive mate- 
rials, the degree and uniformity of heat treatment can vary dramatically with the 
geometric size of the product. 

As a result of these considerations, it is reasonable to expect a close 
quantitative correlation of fracture behavior where details of alloy, process, 
size, and temperature are closely aligned. Where any one of these factors is 
allowed to vary, anomalous fracture behavior can be expected. Although the dif- 
ferences may be rationalized in a qualitative manner, they cannot be assessed 
with much quantitative satisfaction. 

It appears that further insight to fracture behavior is still dependent on 
a continuing compilation of fracture data, until a broad enough reservoir of data 
is available to enlarge the analysis. 

Characterization of Fracture Behavior by 
Stress-Intensity-Factor Concepts 

The severity of the crack-tip elastic -stress field can be defined by the 
stress -intensity-factor concepts of linear elastic fracture mechanics. The 
general analytical formulation of the stress -intensity factor is 

K = S/7 f(a,c,w) , (70) 

where f(a,c,w) is a geometric scaling function dependent on crack and specimen 
geometry. For the specimen configurations considered in this program, the speci- 
fic formulations are, for the compact-type (CT) specimen, 

K = (P/Tw)/r [29. 6-185. 5(a/w) + 655.7(a/w)^ 

-1017. 0(a/w)® + 638. 9(a/w)^] ; (71) 

for the center -cracked (CC) tension panel, 

K = S/F [TTsec(TTc/w)]^ ; (72) 

and for the part-through crack (PTC) or surface-flaw specimen, 

K = S/a [1.21 tt/Q]^ , (73) 

where 

Q = [E(a/c)]2 - 0.212 (S/TYS)^ , (74) 

and E(a/c) is the complete elliptic integral of the second kind. These are the 
basic formulations which were used for evaluating fracture data in this program. 
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As will be seen in the ensuing discussion, fracture instability is not a 
discrete event for most structural materials. Thus, the characterization of dif- 
ferent stages of the fracture process will be identified by different subscripts 
on the stress and crack dimensions contained in the above expressions. 

Crack Behavior Associated With Fracture 

As a basis for comparing of the fracture behavior of various specimen and 
crack configurations, a brief description of the crack extension associated with 
the fracture process is presented in the following discussion. As will be seen, 
there are a number of important bench marks associated with crack extension prior 
to fracture. Any or all of these may be noted in a particular specimen and crack 
configuration. In order to make a rational correlation of the characteristic 

fracture parameters which are derived, it is necessary to relate the important 
bench marks in a comparable and equivalent manner. 

Crack Extension and Specimen Response . - Under a rising load, a fatigue - 
precracked fracture specimen deforms initially in the linear and elastic manner 
shown in figure 23. During this initial stage of loading, the crack extension 
and plasticity associated with specimen deformation are nonexistent or, at 
least, negligible. At some point of loading, a nonlinearity in the specimen 
load -deflection curve is noted and may be attributed to a combination of crack 
extension and plasticity. The degree to which each process prevails could be 
characterized by unloading and marking these specimens; however, this is usually 
not done in the general characterization of fracture. It is only important to 
note that the two processes can and do interact to develop the nonlinearity. 
Finally, after sufficient loading and crack extension, a strain or energy insta- 
bility will develop to fracture the specimen. 

Parameters of Fracture Characterization . — From the previous descriptions, 
there are at least three bench marks to which fracture characterization parameters 
can be referenced. These are 

• The onset of crack extension 

• Apparent fracture instability 

• Critical fracture instability. 

These points are indicated on figure 23 as the points Q, A, and C, respectively. 
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Figure 23, Relation of crack extension and specimen response during fracture test. 



At this point, it is appropriate to point out and repeat the distinction be- 
tween the general concepts of fracture toughness and residual strength. If frac- 
ture occurs in an abrupt fashion after only minimal amounts of crack extension, 
Point C will occur very close to Points Q and A, such that the fracture event is 
essentially a discrete and unique point. The resultant characterization param- 
eter would be a "fracture toughness" value for the material. However, if Point C 
is removed from Points Q and A, all points, as well as the intervening curve, are 
important descriptors of the fracture behavior. This latter behavior is gener- 
ally referred to as residual strength for which one partial, but incomplete 
measure is the "apparent" fracture toughness at Point A, 

The fracture parameters which are associated with these points for various 
specimens are indicated in table 14. 

TABLE 14. 

FRACTURE PARAMETERS ASSOCIATED WITH BENCHMARKS 
OF STRESS AND CRACK EXTENSION 


Stress State 

Plane -Strain 

— 
Plane-Stress and Transitional 

Specimen Type 

CT 

CC 

PTC 

Onset of Crack Extension 

K 

K , K_ ^ 

K , K ^ 


Ic 

q’ Ic 

q’ Ic 

Apparent Fracture Instability 

K 

K 

K 


app 

app 

app 

Critical Fracture Instability^ 

— 

K 

c 

— 


3 . 

Provided that plane-strain criteria are satisfied. 

Generally not monitored in compact specimen and part-through crack specimen 
tests . 


Data Evaluation 

As a characterization of terminal crack behavior, the compilation of frac- 
ture toughness and residual strength data on the subject materials has been 
limited to those stress state and specimen configurations which are most rele- 
vant to the fatigue -crack-propagation studies. These are the quasi-plane-strain 
fracture toughness as determined by compact specimens in accordance with ASTM 
Designation: E-399-72 (ref. 3), quasi-plane stress and transitional fracture 
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toughness, as determined by center -cracked tension panels, and part-through crack 
(or surface flaw) fracture toughness, which simulates natural crack conditions. 
While each of these has a distinct role in the analysis of damage tolerant struc- 
tures, their interfaces are not completely clear because of the interdependent 
complexities of geometric configuration, stress state and basic material proper- 
ties . 

The evaluation of experimental data of any nature involves two basic steps, 
namely screening and analysis. The data are screened in order to assure satis- 
faction of the basic criteria on the characteristics being evaluated. The analy- 
sis, of course, is concerned with computation of the characteristic parameters. 

Screening Criteria . — Within the concepts of linear elastic fracture mechan- 
ics, there are two basic constraints imposed on fracture data to assure the 
characterization of elastic fracture instability. These are frequently referred 
to as the net section stress criterion and the size requirement. The former 
assures that the stress on the gross structural section is dominantly elastic at 
failure; the latter reflects the degree of local plasticity which may be mani- 
fested adjacent to the crack tip. Together, these constraints determine the 
* 

validity of the test as a representation of elastic fracture. 

Net Section Stress Criterion. ~ The criterion which is imposed on fracture 
toughness and residual strength data to assure elastic fracture conditions has 
evolved from experience and, to a large degree, is approximate for each specimen 
type. 

The net section stress is the nominal stress on the uncracked section deter- 
mined in accordance with elementary concepts of strength of materials. It does 
not include the stress concentrating effect of the notch or crack and is used 
only as a simple measure of the nominal stress conditions on the load bearing 
area of the specimen. The net section stress formulations and ratios are de- 
fined in table 15. For the compact fracture specimen, the net section stress 
includes both a bending and tension stress component due to the load eccentric- 
ity. For the center -cracked and part-through crack specimen, the net section 
stress is simply a tension stress on the uncracked area due to axial loading. 


*>v 
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It is important to recognize that in this context the terms “valid” and 
"invalid” refer to the adequacy of the elastic -stress condition and do not 
question the authenticity of the test per se. 



TABLE 15. 


NET SECTION STRESS CRITERIA FOR VARIOUS 
SPECIMEN-CRACK CONFIGURATIONS 


Specimen Type 

Net Section Stress, 

Net Section Stress Criterion, 


S 

n 

S /TYS 
n 

Compact 

-L. mil 

T • w a ^ 

(i-f) 

0.8 

Center -Cracked 

s 

0.8 




Part-Through Crack 

s 

0.9 


(1_^ a ^ 

'' AT w 7 



Size Requirement. — Although the size requirement is not totally independent 
of the net section stress criterion, its consideration arises from slightly dif- 
ferent concepts. In this context, it is used primarily as a thickness require- 
ment on the compact specimen to assure plane-strain constraint of the plastic 
zone. In accordance with ASTM Designation: E-399-72 (ref. 3), this requirement 

is stated as 

T ^ 2.5 (Kj^/TYS)2 . (75) 

It should be noted that while this criterion is also imposed on the crack length 
within the above standard, the previous net section stress criterion is even 

more restrictive on crack length, such that it need not be included here. 

Data Analysis . — The basic fracture data in the form of specimen and crack 
dimensions, loads and stress levels have been analyzed in accordance with equa- 
tions (71), (72), and (73) subject to the above screening criteria. Specifi- 
cally, the combinations of load or stress levels, and crack size dimensions used 
with these equations for the parameters listed in table 14 are, for the compact 
specimen, 

^Ic"®^^q’^^ (76) 

^app ~ ®(^max’^^ ’ (77) 
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for the center-cracked specimen, 


Kq=g(Sq,Co) (78) 

(79) 

\ = g(Sc»c^) . ( 80 ) 

and for the part-through crack specimen, 

Kq=g(S^,ao) ( 81 ) 

^app * ' (®^) 


The analyses have been performed by digital computer using the program 
listed in Appendix G. The output of such analyses are available as a tabular 
format of basic fracture data and associated fracture parameters. For each 
specimen type, the data are categorized by material alloy, product form, thick- 
ness, grain direction, and buckling restraint of the crack edge (if appropriate 
to the thickness). Within this grouping, tabulations are presented by subcate- 
gories of test temperature and specimen size or width. 

Results 

The tabulations of data which have been compiled and evaluated in accordance 
with the previous procedures are described in the following subsections. Al- 
though the formats have been developed to consolidate the data on a common basis, 
there are variations which reflect the different quantities and measurements in- 
volved in each t 3 rpe of test. 

Compact Specimen . — The compact specimen is used primarily to determine the 
plane-strain fracture toughness of relatively thick materials. A sample tabular 
format for the output of this type of fracture data is presented in figure 24. 
Since, at the present time, the initial fatigue precrack length, the 5 percent 
secant offset load and the maximum load are the principal quantities derived 
from such a test, these quantities are presented along with the specimen dimen- 
sions as basic data. The analysis results are presented as toughness values 
associated with the offset and maximum load calculated in accordance with equa- 
tion (71), using the combination of equations (76) and (77), respectively. The 
effective net section stress ratio and size requirements are presented as valid- 
ity checks. Finally, the data source reference is listed. 
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Figure 24. — Sample of tabular output for compact specimen fracture analysis 
results. 
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Center -Crack Specimen . — In that the center-crack specimen is used primarily 
to determine plane stress and transitional fracture toughness of relatively thin 
materials, more experimental quantities are usually recorded. A more expansive 
tabular format for this specimen type is shown in figure 25. The initial fatigue 
precrack length, the critical crack length as reported by the investigator, the 5 
percent secant offset load, and maximum load are tabulated along with the speci- 
men dimensions as basic fracture data. From these items, the offset, apparent 
and critical toughness are computed in accordance with equation (72) and the 
dimensional combinations of equations (78), (79), and (80), respectively. The 
associated net section stress ratios are also presented as validity checks. 
Finally, the data source reference is noted. 

Within the field of each table, the data are categorized and grouped by test 
temperature and specimen width. Following each grouping, where more then one 
valid value exists, an average value and standard deviation are presented. 

Part-Through Crack Specimen . — The part-through crack or surface flaw speci- 
men is used primarily as a direct representative of naturally occurring cracks in 
structural materials in a wide range of thicknesses. As a result, these speci- 
mens and their data can reflect a full range of stress states. A sample illus- 
tration of the tabular format for these data is presented in figure 26. Because 
this crack shape is generally semielliptical in shape, two dimensions, length 
and depth, are required for its description. The initial precrack size, 5 per- 
cent secant offset stress, and maximum stress are presented along with specimen 
dimensions as basic data. The toughness values are computed in accordance with 
equation (73) and the dimensional combinations of equations (81) and (82). The 
net section stress ratio is presented as a validity check. The shape ratio is 
included as an indication of the ellipticity of the shape. Finally, the data 
source reference is noted. 


CONCLUSIONS 

As a result of this study, it was found that large amounts of fatigue, fa- 
tigue-crack propagation, and fracture data can be consolidated for use in design 
applications. Each of these three areas of material behavior were treated sep- 
arately, using large files of pertinent data that were gathered on 2024 and 7075 
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Figure 25. — Sample of tabular output for center -cracked specimen fracture 
analysis results. 
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aluminum alloys, Ti-6A1-4V alloy, and 300M steel. Fatigue and fatigue -crack- 
propagation analyses were limited to constant-amplitude cycling conditions. 

From studies of fatigue data, it was concluded that 

(1) The Walker equivalent strain parameter can be used to account 
for effects of stress ratio. 

(2) A local stress -strain analysis, which uses a computed value 
and a technique to approximately account for relaxation of mean 
stress, can be used to account for notch effects. 

(3) The inverse hyperbolic tangent function can be employed to model 

fatigue curves in terms of e versus log N. for both unnotched 

eq r 

and notched specimens. 

(4) Using the tanh ^ function, it is possible to compute mean fatigue 
curves and tolerance limit curves for 90 and 99 percent proba- 
bility of survival with 95 percent level of confidence. 

From studies of fatigue -crack-propagation data, it was concluded that 

(1) Crack growth curves can be simply and effectively approximated 
using a five-point, divided-difference scheme. 

(2) The Walker effective stress -intensity formulation can be used 
to account for stress ratio effects. 

(3) The inverse hyperbolic -tangent function can be used to model 
crack-growth rate curves . 

(4) Using the tanh“^ function, mean growth rate curves and 90 and 
99 percent probability two-sided tolerance bands with 95 per- 
cent confidence level can be developed. 

From studies of fracture toughness and residual strength data, it was con- 
cluded that 

(1) Consistent fracture characterization can be achieved with stress- 
intensity-factor concepts within a common categorization of the 
details of alloy, process, size, and temperature. 

(2) Correlation of fracture behavior for various specimen types and 
stress-state conditions must be based on equivalent degrees of 
crack extension. 

(3) A broader characterization of fracture data reflecting the influ- 
ence of thickness effects, processing variable grain direction, 
and specimen configuration requires a continued expansion of the 
data reservoir. 
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CYCLIC STRESS -STRAIN DATA 

The method of fatigue analysis developed in this program required the use 
of both cyclic and monotonic stress-strain curves. Using the data source refer- 
ences of Appendix A and information from MIL-HDBK-5B (ref. 1), it was possible 
to characterize the monotonic stress -strain response for the materials of inter- 
est. However, outside of the data reported by Endo and Morrow (ref. 16), 
Landgraf et al (ref. 17), Smith et al (ref. 18), and Gamble (data source ref. 
90), there was no appropriate information available on the cyclic stress-strain 
response of these same materials. To fill this void of information, a limited 
amount of complementary tests were conducted on 2.29 mm (0.09 in.) thick 2024- 
T3 and 7075-T6 aluminum sheet. Specimens were from the same lot of material 
used in a number of previous experimental programs (data source refs. 1 through 
9). 

All tests were performed using an electrohydrualic test system operated in 
closed-loop strain control at a constant strain rate of 4 x 10“® sec“^ . Experi- 
mental procedures were similar to those reported by Jaske et al (ref. 37 ). 
Loading was axial and special lateral guides were used to prevent buckling. 

These guides were clamped about the specimen with a force light enough to avoid 
significantly influencing loading of the specimen. Strain was measured over a 
1.27 mm (0.500 in.) gage length using a special extensometer with a linear vari- 
able displacement transformer (LVDT) as the transducer. Load was measured by a 
standard load cell in series with the specimen and continuously recorded on a 
time-based chart. Load-strain records were made periodically using an X-Y 
recorder. 

Results of these experiments are summarized in table Bl. For each alloy, 
three incremental step tests (ref. 17) were used to develop continuous monotonic 
and cyclic stress-strain curves up to 0.01 maximum strain (see figs. Bl and B2). 
To see if the cyclic stress-strain curves from the step tests could be used to 
predict cyclic stress-strain response under constant -amplitude strain cycling, 
seven specimens of each alloy were tested under constant -amplitude loading. For 
three tests the strain ratio (algebraic ratio of minimum to maximum strain) was 
equal to -1.0 (i.e., the mean strain was zero). A positive value of mean strain 
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was used in the other four tests - three were with a strain ratio of 0.0 and one 
was at a strain ratio of 0.5. 

In all cases, results from the constant-amplitude tests were close to those 
predicted by the cyclic stress-strain curve from the step tests (figs. B1 and 
B2). Thus, it was concluded that these cyclic stress-strain curves could be used 
to describe the stable stress-strain response of these two materials. 

Unpublished cyclic stress-strain data have been generated on 300M steel and 
annealed Ti-6A1-4V alloy during in-house studies conducted at BCL. Experimental 
procedures were the same as those described earlier, except that a 0.635 mm 
(0.250 in.) diameter, 1.27 mm (0.500 in.) gage length specimen was used. Cyclic 
stress-strain curves for these two alloys are presented in figures B3 and B4. 
Samples of the titanium alloy from the transverse (T) direction and from elec- 
tron-beam (EB) welded plate cyclically hardened. Whereas, samples from the lon- 
gitudinal (L) direction cyclically softened. The cyclic curve shown in figure 
B4 is for the L direction and the monotonic curve was estimated from published 
data (data source ref. 70). To show the wide variation in cyclic stress-strain 
behavior of this alloy, data from Smith et al (ref. 18) are presented in figure 
B5 and data from Gamble (data source ref. 90) are presented in figures B6 and B7. 
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STEP -1.0 0.0206 max -- — 30-37/40 

CA -1.0 0.0201 0.0056 1 050 -- 292 














TABLE Bl. -RESULTS OF CYCLIC STRESS-STRAIN TESTS AT A STRAIN RATE OF 4 x 10“^ sec“^-Concluded 
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Figure Bl. - Cyclic stress-strain behavior of 2024-T3 aluminum 
sheet. 
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Figure B5. Cyclic stress -strain behavior of solution-treated and 
aged (STA) Ti-6A1-4V bar, data from Smith et al (ref. 18). 
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Strain 


Figure B7. — Cyclic stress-strain behavior of annealed T1-6A1-4V 
bar, data from Gamble (data source ref. 90). 
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STEP-BY-STEP APPROACH TO ANALYSIS OF FATIGUE DATA 

The following outline provides an illustration of the procedure which was 
developed for the consolidation and analysis of fatigue data. Each part of the 
analysis is broken down in a step-by-step fashion so that the details of the 
process might be more clearly defined. Possible simplifying assumptions are also 
noted. 

A. Material Data Acquisition 

(1) Constant -amplitude smooth-specimen fatigue data (preferably obtained 
over a range of stress ratios or mean stress values) are required to 
optimize the equivalent strain m factor. If only notched data are 
to be analyzed, an optimum m value must be estimated. 

(2) Constant -amplitude notched -specimen fatigue data make up the second 
part of the data file for a given material. Again, data generated 
at various values are desirable if a reasonable value for p is to 
be determined. 

(3) Monotonic and cyclic stress-strain information for the investigated 
material, heat treatment and product form are required for both the 
smooth and notched fatigue specimen analyses. In order to perform 
the analyses, it may be necessary, in some cases, to develop reason- 
able estimates of the cyclic and monotonic curves on the basis of 
available information. This should be done with considerable cau- 
tion, however, since effectiveness of the entire analysis is based 
on appropriate material property definition. 

B. Smooth Specimen Fatigue Data Analysis 

(1) Calculate maximum nominal stress and alternating strain values for 
the complete data set. 

(a) If the input is in terms of maximum strain and 

strain ratio, calculate alternating strain values 
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for each specimen according to the following 
expression: 


e 

a 


^max 


(1-R) 

2 


(Cl) 


and determine maximum stress values through 
usage of the cyclic stress -strain curve [eq. 

( 10 )]. 

(b) If input is in terms of maximum stress and 

stress ratio, calculate values of alternating 

stress in the same manner as shown in equation 

(Cl), substituting values of S and S for 

a max 

values of e and e . Then determine alterna- 
a max 

ting strain values according to the cyclic 
stress-strain curve [eq. (10)]. 

(2) Convert cycles to failure to log N^. 

(3) Set limits within which the optimum m value is likely to occur. 

In most cases, limits of 0.30 and 0.50 would bound the optimum 
m value. 

(4) For a particular m value, calculate values of equivalent strain 
for the entire smooth data set according to equation (9b). 

(5) Fit a third order polynomial [eq. (29)] to the set of calculated 
equivalent strain values through a least-squares regression process. 
(The polynomial, rather than the tanh“’- function, was used in computa- 
tions involving an optimization on m because the quality of fit for 
the polynomial was related solely to the degree of data consolidation. 
Quality of fit for the tanh“’- function, however, was somewhat dependent 
on the specific function limits which were used, and these limits could 
not be appropriately determined until a specific m value was chosen. ) 

(6) Calculate and record the total sum of squares and the sum of squares 
of deviations for the regressed equation. Then calculate an value 
according to equation (4). 

(7) Repeat steps B4 through B6 for complete range of m values and select 
the best m value in terms of maximum R^ . (This procedure is almost 
necessarily handled by a computer using an iterative process.) 
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(8) Using the optimum m value (or a satisfactory approximation), perform 
a least-squares regression to determine the best fit inverse hyper- 
bolic tangent equation for the investigated data set. If the data 
cover the entire fatigue life range from 10 to 10® cycles, establish 
the range of the function according to equations (36) and (37). 
Otherwise reasonable values must be specified for these limits. 

(9) Apply the weighting function W(e^^), [eq. (42)], to the residuals 
and determine whether the modified residuals are sufficiently uni- 
form throughout the range of data. It may be desirable to examine 
normality through construction of a frequency distribution plot of 
the residuals. This is done by plotting the frequency of occurrence 
versus the magnitude of deviations from the mean curve. 

(10) If step B9 is completed satisfactorily, probability of survival 

lines can be constructed according to equation (7), and the result- 
ant curves may then be displayed graphically along with the individ- 
ual equivalent strain data values. 

C. Notched Specimen Fatigue Data Analysis 

(1) As in the smooth specimen fatigue analysis, nominal values of maxi- 
mum stress and alternating strain must be calculated in the analysis 
of notched specimen data. Steps Bla and Bib are applicable. The 
cyclic stress-strain function [eq. (10)] is used in both cases. 

(2) To calculate local mean stress values according to equations (23) 
and (27a), it is also necessary to calculate monotonically based 
values of nominal maximum and alternating stress and strain. The 
same stress-strain equation [eq. (10)] is used for these calcula- 
tions as was used in step Cl, but monotonic parameters (as in table 
6) are required. 

(3) Convert cycles to failure to log N^. 

(4) Set limits within which the optimum p value is likely to occur. 

(In most cases, limits of 0.00 and 0.03 would bound the optimum 
p value.) 
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(5) For a particular p value, calculate values for the fatigue con- 
centration factor according to equation (20). Knowledge of the 
notch root radius is required for this calculation. 

(6) Calculate an estimate of local alternating strain values accord- 
ing to the following expression: 


= K-e 
f a 


(7) Then calculate approximate values for local maximum stress using 
the following relationship: 


a = a + cr 
max a m 


(C3) 


It is not possible to calculate directly because cyclic 

plasticity allows mean stress relaxation that decreases stable 
local mean stress values. Equation (C3) is an approximate means 
of accounting for that relaxation for constant -amplitude loading. 

(a) In this equation, a is calculated from equation (10) 
using values of local alternating strain determined in 
equation (C2). 

(b) Values for a in equation (C3) are found according to 

m 

equations (23) and (27a). 

(i) In equation (27a), the value, e' /2, is equivalent 

max 

to a monotonic value of local alternating strain. 
This quantity can be determined by multiplying the 
monotonic nominal value of alternating strain, 
found in step C2, by the fatigue concentration 
factor which was found in step C5. 

(ii) Similarly, the magnitude of the monotonic local 
maximum strain used in equation (23) can be deter- 
mined by multiplying the mono tonic maximum nominal 
strain found in step C2, by the fatigue concentra- 
tion factor K^. 

(8) For a particular p value and for an optimum m value found in part 

B, calculate values of equivalent strain for the entire notched 

data set according to equation (9b). Values of a and e calcu- 

max a 

in steps C6 and C7 are used, however, in place of S and e . 

max a 
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(9) Fit a third order polynomial [eq. (29)] to the set of equivalent 
strain values in the same manner as in step B5. 

(10) As in step B6, determine TSS, SSD, and R® values for the regressed 
equation. 

(11) Repeat steps C5 through CIO for the complete range of p values 
and select an optimum p on the basis of a maximum R^ value. 

(Again, the computer is almost essential for this operation.) 

(12) Using optimum values for m and p, determine a best-fit inverse 
hyperbolic tangent equation in the same manner as in step B8. 

(13) As in step B9, apply the weighting function, examine the residuals; 
and if necessary, construct a frequency distribution plot of the 
residuals . 

(14) If step C13 is completed satisfactorily, construct probability of 
survival lines as in step BIO. 
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RESULTS OF CONSTANT -AMPLITUDE FATIGUE DATA CONSOLIDATION 


The collection of figures in this appendix displays results of the constant- 
amplitude fatigue data consolidation effort. One figure is presented on each page, 
and each figure consists of two related plots. The upper plot shows the consoli- 
dated data along with the regressed mean curve. Below the mean curve, are the 
calculated 90 and 99 percent statistical tolerance curves, respectively, which were 
established at a 95 percent confidence level. The lower plot illustrates the 
pattern of weighted residuals for the consolidated data shown in the upper plot. 

As explained earlier (see page 30), the abscissa is in terms of actual fatigue life. 
The residual plots were included to provide a visual indication of whether the 
statistical requirements of randomness, zero mean deviations, and uniformity of 
variance were satisfied for the consolidated data. 

Figures D1 through D4 represent consolidated data for 2024-T3 and 2024-T4 
aluminum. Figures D5 through D9 consist of consolidated data on 7075-T6 and 7075- 
T651 aluminum. For both series of aluminum, the data consolidation was substantial. 
The best collapse of data (R® = 99.6 percent) was obtained for the unnotched 7075- 
T6 clad sheet material, where all data came from a single source. In all but one 
case (fig. D8), consisting of both notched and unnotched aluminum data, the R® 
value was 94 percent or greater. Other than the noted exception, unnotched data 
were consolidated better than notched data. Nonuniformity of variance was of some 
concern in several cases (figs. D2 and D6), but this problem was not due to 
inadequacy of the weighting function; it was due to layering of data from different 
sources in the high cycle fatigue range. This layering made it impossible to 
account for data in this region as effectively as data in the lower cycle regions 
where no such layering was evident. The nonuniformity of variance was not con- 
sidered to be severe enough, however, to make the construction of tolerance limits 
inappropriate. 

Results for the 300M steel fatigue analysis are presented in figures DIO and 
Dll, The R® values for both curves were not as high as the values determined for 
the aluminum alloys, but the overall data collapse was considered good since the 
inherent data scatter for the two data sets was quite large. 

The T1-6A1-4V alloy data, displayed in figures D12 through D23, were the most 
difficult to analyze and provided the poorest results. The difficulties were due 
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to two major factors. First, the titaniutn data file consisted of a large number 
of different product forms and heat treatments, and although an attempt was made 
to develop accurate monotonic and cyclic stress-strain data for each variation, 
only a rough approximation of these curves was possible in some cases. Secondly, 
the inherent scatter in most of the titanium data was very great, making a 
consolidation effort difficult. Despite these problems, R® values exceeding 80 
percent were obtained in figures D12, D17, D19, D20, D22, and D23. The best 
results were found for the T1-6A1-4V data in the solution-treated and aged 
condition. For cases where the data consolidations were not acceptable as shown 
in figures D13, D14, D15, D16, D18, and D21, the values of R® were below 80 per- 
cent. Such curves cannot be used for design applications and are included in 
this report only to show how poorly the analytical procedures worked in some 
instances. Until more experimental information is obtained upon the cyclic 
stress-strain and fatigue behavior of Ti-6A1-4V alloy, it will not be possible 
to refine the analytical procedures to account for such anomalies. 
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b. Distribution of weighted residuals. 


Figure D3. — 2024 -T4 Bar and rod (unnotched). 
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Figure D14. — Annealed Ti-6A1-4V bar, extrusion, and forging [140 ksi TYS 
(unnotched )] . 
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a. Consolidated fatigue data with mean curve 
and 90 and 99 percent survival lines. 



b. Distribution of weighted residuals. 


Figure D15. - Annealed Ti-6A1-4V casting (unnotched). 
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Figure D19. ~ Annealed Ti-6A1-4V casting (notched) 
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a. Consolidated fatigue data with mean curve 
and 90 and 99 percent survival lines. 
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Figure D23. — STA Ti-6A1-4V casting and plate (notched). 
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APPLICATION OF FATIGUE -CRACK-PROPAGATION RATE ANALYSIS 

An actual structural element functioning in a real world environment is 
usually subjected to an extremely diverse loading spectrum. Large fluctuations 
in maximum stress, stress ratio, and frequency, not to mention temperature and 
environment, are likely to occur. 

A crack extending in a structural element in such a loading environment 
will propagate in a complex manner. It is this highly complex situation that 
the designer must consider when structural life predictions are made. Clearly, 
some simplifications of the actual loading situation must be made for the pur- 
poses of analysis. Crack-growth-rate analysis is often used to predict life. 
Such an approach provides a model of the damage process by which cracks grow 
under constant -amplitude cyclic loading. Once the characteristic model for a 
material is obtained, it can be integrated to yield life predictions. Applica- 
tion of the rate analysis approach developed in the program is discussed brief- 
ly in this appendix. 


Fatigue -Crack-Growth Rate Model 

Previous discussions in the main body of this report have examined how 
crack-growth data can be converted to a crack-growth-rate format and subsequent- 
ly represented by an analytical model of fatigue-crack propagation. Using this 
approach, data from a large number of sources, covering a wide range of parame- 
ters, can be combined to characterize a material quite thoroughly. As a result 
of this investigation, the inverse hyperbolic tangent model has been proposed 
as a suitable tool for data analyses. 

Rate Integration 

Once a satisfactory composition for the governing differential equation 
has been developed, the second part of the fatigue -crack-propagation problem, 
integration, becomes important. Some of the possible integration procedures 
are presented here in schematic form. However, implementation of these proce- 
dures was outside the scope of the present work. The general expression, 

da/dN = f[K(a)] , (El) 
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is a first order, linear differential equation. Although the fatigue-crack- 
propagation law is in the form of a rate equation, it is essential that the 
designer be able to calculate crack length for the case where a given number 
of cycles have been applied under specified loading conditions. This calcula- 
tion can be accomplished by using the result of integrating equation (El), 


da 

f[K(a)] 


N, 


= f* 


dN 


(E2) 


N, 


Mathematically speaking, this integration is representative of a class of solu- 
tions of differential equations known as "initial value problems". In other 
words, it is desired to find the final crack length after a certain number of 
cycles have elapsed, starting from initial values for crack length and cycle 
number . 

1\70 general methods of solution of the initial value problem are available 
The first, and most straightforward method consists of carrying out the integra 
tion indicated in equation (E2) to yield a closed-form solution. The closed- 
form solution has the advantage of being a concise equation from which the num- 
ber of cycles required for a crack to grow to a given length is easily computed 
Unfortunately, it is not always easy, or even possible, to perform the integra- 
tion of the differential equation. The second method of solution involves the 
use of a numerical integration scheme such as the Runge-Kutta method. 

Hoskin (ref. 30) gives some closed-form soltuions to the most common 
fatigue -crack-propagation models where the variables can be readily separated 
and integrated. Consider, for example, Paris' Law, 

da/dN = C(AK)'' . (E3) 


Assuming K = 
max 

negligible, the 


S /rra for the case of center-cracks where width effects are 
max 

integration for n' > 2 is given as 


N = 


r. 


max 



n^(l 


->T(¥ 




(E4) 
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A similar procedure may be carried out for the Forman -model differential equa- 
tion. 


da/dN = 


C(AK) 


(1-R)(K -K ) 
c max 


(E5) 


yields for n^ > 3 


N = 


l^max 


""'“V (n'-3)/2n 


:®c/^i) 

7n^ 


n ' _i I 

T 


l-(a^/a) 


L"i 


1 

■ (n^-3) 


^ (n" - 3/2) 
^i 


(E6) 


When the K expression involves more complex algebraic or transcendental 

functions for which a closed form solution often cannot be achieved. Such is 

the case when width effects are not negligible. For example, when a width cor 

rection factor term is used, the expression for K becomes 

’ max 

K = S y TTa sec(— ) , (E7 

max max^ \ w / 


and the integration is much more complex. 

The alternative to closed-form solution of equation (E2) is numerical inte- 
gration, which necessitates the use of a digital computer. Engle (ref. 38) for 
example, utilized the Runge-Kutta numerical integration scheme in Program CRACKS 
(which, incidentally, accommodates variable -amplitude loads as well as the more 
elementary constant -amplitude cases). This is typical of the type of numerical 
solutions used today. Since the inverse hyperbolic tangent model has been based 

on the above form for K , equation (E7), the integration process becomes exceed 

max 

ingly complex. A closed-form integration of the model does not seem to be practi 
cal. However, the model may be integrated numerically. The expression may be 
put in the form. 


Ng-Nj^ = 10" = 


Cain 


log 




K K 
o c 


(l-R)^*” S^ TTa sec rra/w 


da 


log . log(^ 


(1-R^ TTa sec Tia/w 


K K 
.....q.C- 


152 



APPENDIX E 


which is suitable for numerical integration. Several integration schemes are 
directly applicable to this expression. These include iterated Gaussian quad- 
ratures as well as the Runge-Kutta method mentioned earlier. 

Summary of the Life Prediction Procedure 

The procedure for predicting the life of a structure, neglecting variable- 
amplitude loading effects is outlined below, 

(1) Crack growth data on the material of interest is collected and 
combined. 

(2) Crack growth data are converted to crack -growth -rate values by 

means of the five -point divided difference scheme. Values of 

K ^ are calculated, 
max 

(3) The inverse hyperbolic tangent function is fitted to the data 
to yield regression and optimization coefficients, as well as 
statistical parameters. 

(4) Initial and final crack lengths, the stress levels, and mean 
stresses for the structures to be investigated are specified. 

This information may be presented in the form of load spectrum. 

(5) Finally, the values specified in item (4) are used to integrate 
the characteristic fatigue-crack-propagation model for the 
number of cycles required to extend the crack to the specified 
final crack length. 
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RESULTS OF FATIGUE -CRACK -PROPAGATION ANALYSIS FOR FIVE MATERIALS 

Figures FI through FIO present the results of analyses performed on 2024-T3, 
7075-T6, and 7075-T7351 aluminum alloys, 300M steel, and Ti-6A1-4V alloy. The 
composition of these data sets was previously described in the report. (See 
pages 79 and 80. ) Each material is characterized by a crack-growth-rate curve 
and an accompanying plot of the distribution of residuals as a function of actual 
crack propagation ate. 

The crack-growth-rate curves consist of the experimental data plotted on a 
= (l-R)”*K ,nav basis. The best-fit regression curve is represented by the 
solid central line through the data points. On either side of the mean curve 
are the 90 and 99 percent tolerance limits, with 95 percent level of confidence. 
Coefficients of the inverse hyperbolic tangent model and the tolerance limit 
formula were presented previously in Table 12. 
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Fatigue Crack Propagation Rate, d(2c)/dN, m/cycle 
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Stress Intensity Factor Range, Kgff, 

Figure F3 Fatigue-crack-propagation-rate curve for 7075-T6 
aluminum alloy. 
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Fatigue Crack Propagation Rate. d(2c)/dN, m/cycle 




Fatigue Crack Propagation Rate, d(2c)/dN, m/cycle 
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Figure F9. - Fatigue-crack-propagation-rate curve for Ti-6A1-4V alloy. 
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COMPUTER PROGRAM FOR FRACTURE ANALYSIS TABULATION 

A listing of the computer program, FRACTAB, for the analysis of fracture 
data is presented in this appendix. The data input formats for this computer 
routine have been presented in the subsection on Data Recording and Storage. 
The typical printout formats have been illustrated in figures 24, 25, and 26. 
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C 

C 

c 

c 


c 


Tow ^NGUTSh PLACt S^lTC'i»2* CA«D IN S3U^r>i i)£CA FOR 

LNG.ISH OJTPlIlt O^MlT SWlrCHt2. CARD FOR METRIC OJirPUT 

_FOR METRIC PuaCE SWITCH* 2. CARD IN bOURctfOECK FOR 

METRIC OUTPUT* OMIT S'^irCm2, CARO FOR ENGLISH OUTPUT 


pr»oooi_ 


,.Of>OOP3„ 


PROGRAM FRACTAL ( I NPUT,QJTPJTfTAPE5=INPUT,TA^£6aJuTPUT) 

_ DIMENSION KOdEI («>:)) I KUDEp(dO)* KOMENT(B)* SUMO), SS30), Xn(3)i 
lAVGO), UEV(3)* X<VAP(3)# S’LAGOJt RATI0(3>* ^ii?) 


" lo 
20 
30 
AO 
50 
60 
70 
BO 
90 
100 
llO 


PJ0005 
__ 5 i *0006 


00 10‘ NslfSl 

_HEA3 (5»455> <OUEi(N) ,iOO"?(N) 

030O1S 

POOO^O 

. lOV- , 

IF t EoF «5) 1 03f 1 C 0 
CONTINUE 

00002? 

000024 

\M- 

5T0» 1 

CAUu' ss«rcH n.is<»i) 

rnoojb 
000030 


CAL. SSwiTCH (gfiMET) 
_PRl>^T 4t>0 

000034 
-__0 00035. 


L1M=N 

_THICXar-. 


C 

- c 

NEAOJ^NO STORE NE<< TITLE* E-S£. STOP RON 

000037 

OijOOSS 

Uv> 

REA3 IS*'*bb) IJN 1 T * I T A * I T^, I TC* THICX 
IF (EoF »'■}) 305* 1 1 5 

036060 
060066 

lib 

HEAD (5.S30 ) foment 
READ ( 5*540) rniC<*lTYPE 

000076 




IF UhlCX.EQ.O.) THIcX = THtC<i 
O u Tn 325 


C 

MOVE INPUT RECORD BACK ONr: CARO 

066101 

120 

c 

OACKSPACE 5 

65610, 

c 

1?5 

'XERO SUMS FOR NEXT data 5 ISSET ' 

DO 136 N=_li3.. _ _ _ . 

rujQios 
_ rhjQlOfe 


SUM(N)=5, 

SSO(n)=u« 

000107 

TiV~ 

c 

XN(N)=r, 

666112 

c 

135 

read NE^r DATA CA=1D 

READ (5*545) lOE.NT * ThI C< . I OTH* CSUBO* SS J 82* SS0B3, Y0CCL*PRCCL,STMAX 

66015^ 


l*TY 5*fuS*5PECL*lTYPE.ISPt:,lfEMP»IREF' ' ~ 

IF (EoF*5) 275*140 _ 


c 

TEST FOR BLANX CA?0 (WIOTm» 0) 

nooibi 

l4u 

c 

If (wIOfH.EO.o.) 30 TO 2/6 

6661 62 

c 

ff'SwifcHi Is ON* replace ref c~ooe 

IF ( JSwl .EO.?) SO TO 155 

666164' 

6661 66 


Do 14S N=1 *LI 4 

IF (IHEF .EO.KOOEI (N) ) SO fO 150 

6661 7o' 
6601 7? 

lAb”” 

'continue 

SO TO 155 


000173 15u lREPsKODE2(N) 


A 
A 
A 
A 
A 

A .120 
A 130 
lAO 
150 
100 
170 
180 
l90 
200 
210 
220 
230 

2H0 

250 
200 
270 
2b0 
290 
300 
310 
320 
330 
3A0 
300 
300 
370 
380 
390 
400 
AlO 
420 
A30 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
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aoQiis 




6661 75 ~ 

c 

c 

1 55 

CHECK FOR MISSING INPJT DaTA 
IF (TYS.LT.I.) go to 390 

— ■ 

DQQaftJ 

C 

(3LANK OUT FLASr'Kf' ♦ RATIj, 
_yQ J.tl.Nsl »3.. 


at} 02 o? 

^666203 


AKVAR (Ml aj. 
«ATIO(N)=g, 


6iiC204 

I 6 u 

_.C 

KFLAU(N1=1H 

CHECK INPUT aNO type of OjT’OT 
IF (MNIT.EU.l) GO To 165 
_1F (ImFT.lO. 2) CO TO l 8 ii 

OESIREO ANO '<AKfc aN.Y CHANGES 

606210 
606212 


000214 
»>20214 __ 

165 

60 TO 1 7u 

..IF (ImfT.EQ.?) CO TO 175 


666216 

C 

T30 TO Ido 

.Switch MclRIr. JNITS TO ENGLISH 

.JNIJS 

066217 
— 666222 

170 

T HI CKr Trl 1 CK/n . 0254/ 1 nou 
W 1 0! H = w 1 UTH/p ♦ 0 254 / i 0 00 

0^6225 


CSJ:3rnCSLlriu/n, 025<» / Inoo 



_p3Q23o . 
f>3023l 
0 j0235 
fi j0237 
.i002<tfL 
riHo2«>i 

_Oro2‘i4 
6 o02‘»7 
. 000250 
000251 
_go0253_ 

0:16260 

000261 

_55o264 _ 

or.0267 
6.102 I? 


C 

175 


0J0273 

0002^7 . 

000301 
?J?3U2 

5.56303 
66630^ 

00031 1 

656312 . 

000313 

220317_ 

0. *16 323 


.,Au=5.‘1‘V5 

IF ( ITYPi:.EU. i )" A 3 = 4. 43 " 

SSU32 = SSU:i2/AD 
SSU33 = SSU"t3/A3 
_STM4;^ = STHhx/a3 _ 

VUC:l=V 0CCL/6. 0254/1 000 
FHC:L = P-^CCL/p.u254/l,ioo 
lrS=TYS/b.895 
1 Jb=rUS/6.895 
ITE'^prl lE^lpii^./D.^Sa. 

_St'c:L = SPHCL./6. 0254 / 1QOO_ 

(jO to Ibo 

SWITCH English units to MrT^ic jnits 
THI:k=TH1CK»o. 0254*1000 
w1lMH=wIUTH*o.o254*i6oo . 

CSU3f=CSU '10*0.0254*1000 
A0=5 • H95 


IF (ITyPE.EO.I) AD=4,45 
SSU32=S5JH2*A0 _ 

SSUb3 = SSO'l3*AD ' 

_ STM4X = STMAX*aD _ 

V 0 CCi. = V0CCL*0 . 0254*1660 

_PPCCl = PHCCL*:', 0254 M 600 

TYS=TYS*0.H9S 
, 1 JS= ruS*t).895 
I TEwPs ( IIEMP-32* ) ‘‘5.’/9. ' 
SPE:l = SPECL** 6 . 0254*100 0 
Co lo 150 


C 5RAVCH TO APPROPRIATE tOJATIO'!' 

c 1 __#compact tension* _ _ __ 

C 2 *CEMTER CRACiX* 

_C 3 _.*S'JrFACE F.An - 

C 4 itSURFACt F.AW - 

C 5 *OOU tjL£ Ca N T1 Li V E R JJJ a w * 


6 *NOTCH BtNO* 


570 
5H0 
590 
600 
610 
520 
63U 
6<»0 
65q 
66(1 
670 
.. 6M0 
69o 
700 
7l0 
720 
730 
740 
750 
760 
77o 
7Bo 
790 

900 

Blu 
B2o 
83o 
940 
350 
860. 
8 To 
8130 
890 
900 
9i0 
920 
930 
940 
9bo 
960 
970 
960 
990 
*1000 
*10lO 
AIO2O 
A1030 
A1 040 
A1050 
A1 060 
A1070 
AlOBO 
A1090 
AllPO 
AlllO 
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nr>p3^3 

cjoJ24 

B25326 

0.10330 

000 33? 

In0333 

-6oo3'*2. 


00C354 

-5563&‘>_ 

0003bft 

_0i0357_. 

- 000 ^ 00 - 

0i)04?<f 


C ScLiCT CRITICAL ST.^ESS •'\M3 
l.ao . SCHIT = SSUH3 
CChlTsPHCCL 

IF (ScNi r.Efj.u.) scmi=ssi3i. 

IF (CCPIT.EQ.O.) CCHTTsv'UcC.I 
IF UTY pE. E (J . SL-O a JQ_3^» 


IF ((nrPE*5)/b.NE.n STOP 
GO TO ( 1»S«2?()<243*3 ‘»o* 3*^0* 3.9o> ♦_ ITYPE... 


_c Compact tenston 

165 IF (CSUHo.fcU.O.) 30 TO 3^* 

1F_ ( TMiCK .to , n .) 5 0 TO 39* 

AsCSUHr/wluf H 

XKVAPH )»SSUrt3/wIDTH/TMlC<»SQ!U(CSJd<J>*(29.6-165,5*X*655.7*X*X-l6l 


17.*X**3*t>36.9»A*»A) 

XKVAS(3)sSIMAX/wi0TH/THlC<'»Sai1T(CSJBn)*<29.6-i65,5*X*655.7*X*X-lnl 

17.*X»»3+h3M,9»X*»4) 


-XU LL=i2.3" (X<y AH U >y.T!f S).tt<»2J- 


XH2) = (2.3«(x<var(2)/TYS)»»2) 

030^26 IF ( IHtT .Nt.p) GO T0 _j95 

000*3o DO ) V* N=i »2 

_000432 XKVAR (N) = XKVAH(N)*lonO*50^T(O.OOl) 

0..0437 IHo XI ( Ml = (2.5* ( x<vAr ( N) /TVS) »»2) »1 000 


. OCOAA6 

p56A5o 


C 

195, 


CHECK 2.5(K/TyS)«*2 r,E 6 ”|H CSJ60 
_RATlOn )=X/1 hICk 

RAT10(2)=X/C5J6a 


- A1120 
A1130 

- AllAO 
A1150 

. AUOO 
A1170 

^AllSO 

A1190 

*1200 

A1210 

A1220 

A1230 

_..*12A0. 

A1250 

A126U 

A1270 

*1200 

A1290 

_.A1300. 

A1310 

- *l32o 
*1330 
AI34U 
A13S0 
A1360. 
A1 370 
A1300 
*139o 
*1 AOO 
AlAlO 


OO0A73 

ooosoo 
goo50A_ 
000507 
,936^1 3_ 
090520 ” 

(1.10523 


OwOB23 

J;0525_ 

900526 

.0.3o527_ 

000532 

_n 0.9 5 34 . 

000536 

000540 


c 

-C _ _ 

205 


IF (SSUH3.GT.0.) HATI0(3)=SI'1AX/SSU83 

DO 215 N=),2 

IF (XI (IM) .GT. THICK) GO TO 205 
.IF (XI (fj) .GT.CSOBO) 30 To 2i)5 _ 

IF (iM,FiJ.2) SMET = 5MET*STMAX2SSj63 

JF (SMFT.3T.TY5*.5) 5O_TO_20b 

IF (HATi()('3') .Gl.l.i) GO To. 20'S 

Go To 210 _ 

IMVALIO 

kflag(ni = i’h* ’ 

.H0ASI=2_ 


C CaLCOLATE net 5ECTI0M STHrSsES 4NI) COMPARE XITH YIEUD 

P-0 P 45 1 I FL J 55 U H 3 . G T . . i . _S N E T = 5 S U 3. 3 7 T 4 1 CK/(WI OTM -C S JB 0 ) »(1 .*3 . * U I OJ H T C S U B_. A 1 4 2 0 

... 1 0 ) / ( * I D TH-C5U4U ) ) AI43O 

0A0465_ IF ( I;iET. \j£.21 Gq Tq 2oO . _ _ A1440 

0.10467 Sn£T=SN£T*1 000 A1450 

0-1C471 XX = Smft/TYS __ _ A1460 

C Aj 47 (j 

COMPARE Max load XITh 5 P-HCEMT OFFSET LOAD- SHOULD NO T EXCEFD 1.1 A14BU 




215 

_C 


GO To 215 
XN( M) sXN (N) ♦i’. 

SUM(N) sSUM (Nl ♦XKVAH ( nT 
SSU(N) =SS'J (N) ♦XKVAR (N>.**2. 
COnT InUE 


A1490 
. _ . *1500 
AI 5 IO 
A1520 
A1530 

A1540 

A1550 
A1S&0 
A 1570 
„_*15B0 
A1590 

A16OO 

A1610 

*1620 

A1630 

A1640 

A1650 

Al()60 

A1670 
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000541 

t^.obi i 
600613 . 


w«nE (b»^>00) 
(fi) iKFuA 3 
Ulf‘jib = Ll;4Li>i 
_ .bO T0..13S . . 


c 

c 

t£.f:lJEii_ crack 

000613 

220 

Xsi .57'rt*CSUR6 

boiifelS- 


... IPCI=8* . . ., 


C 

COMPUTE STRESS 


PKni-DNC-LIWl-UF- INPUT/UjT?ur-..JATA 


0A0616 


0l»062A 


5 jC 63 ? , 

0 .j 6 <J 40 

05 Cb 52 _ 

0 ') 06 b 4 
000666 

0 jo 700 „ 

500 ^ 0 ? 
OijC 7 . 04 .^ 


1 *APPAKCNT>s 

IF (CSUBO.QT.5.) PATl0(e!)3Sr'iAX/TYS/(l.-CSU30/«ll0TH> 

2 ~#CPIt1CALx ■■■ 

JF (CCWIT.GT.O.) PATlO(3)sSTMAX/TYS/(I.-CCHIT/Wlr)TH) 

3 xoFFstrjt 

_IF (CSUBU.i^r.O. ) PATtO(1)sS:RIT/TYS/(1.-CSU3«/i»IoTHL 
iF (X.GT.i).) XKVAP (2) =SrM4X»S3-U (X/CuS(X/Wl JIM) ) 

Jf _ U . 0 T vO . ) _ X Ay A i ( I ) j: 5 C Hi T itS 3HT(X/C0S < x H > i 


x=) .57?8*CCHI r 
IF ( X . G T . rt . ) X K V A H ( 3 ) = 6 V Mft X »S3H T ( X_/COS ( X/w I DT M )J_ 


226 


000711 

D' 307 l 3 

^8;»0716_ 


C 

230 


IF (IMET.NE.P) GO T0.230 
00 226 N=) i3 

XiXvAP(N)=XKVAH(N) »SQPT (0.^01) 
CHECK M-S stress PATTo i-F .3 

00 24* N=1 »3 ^ _ 

IF (XKYttH(N).U(.l.) GO TO 2<»0 
IF (RATIO (N) _G0 To 235 


A] 6UU 
AlS'^O 
A1700 
A1710 
A1720 
A1730 
..A17^0 . 
A1750 
*1760 
41770 
A17B0 
41730 
A1900 
41810 
41820 
AiB 3 o 
AI 8 aO 
AiebO 
A18O0 . 
4 i 87 o 
4i880 
A1830 
41900, 
A1910 
_A1920 
41930 
41940 
A1950 
A196u 


C 41970 

C INV4LID 41980 

000720 KFLAG(N)=1H« 41930 

000 722 N(7A5T = i _ . . A^qOO 

000723 GO fO 240 A2QI0 

C VALID ~ - ^2030 

pOQ72 3 235 X(4 ( M) =XinI ( ■'I) ♦ i , A2Q40 

000726 SUM (D ) =SU -1 (tvl) *XNV AH ( N) A^QbU 

c5673o . SSU (rg) =SSJ (N) *XKVAR t^) _ 42 q60 

noo73? 24 o cuTgiiNUE 

_ JC __ _ 42080 

C PKlvT ONE LINE OF INPUt/OiTPUT 34TA 42030 

530734 1T = 2 42100 

000736 HRI(£ (6,405) 1DENT,ITF.NP,TYS»TJS, THICK, (ilDTH,CSj30»CCRlt,SCRit, ST 421l0 

... IHAX, ( (XKVARdi) ,KFi.* 8 (N) ) , v(= 1 , 3 ) * HAT 1 0.< I > fHAT.IOlgJ ♦HaTI»( 3) fiREFi . 4 ;>l 20 

noibus lim- s=Lii,Es*i a2i3o 

001007. GO TO 135 A2l4u 

C 42150 

e SlJRFACe FiAW 4216(J 

001007 245 IF (SOECL.E( 7 . 0 ) GD To 330 A2l7o 

_ Poi&lO IF (CS(J60.E(3.0) GD TO 330 _ _ _ _ . 4j>18U 

oolbli IF (STMAX.Eq.O) Go To 330 42130 

001012 1 PCT = 9^ . . __ __ 42200 

C A2210 

,, £ !;iatL=LUI£-4SgEci_ a ND ST R es s ratio s .4 222 0 

001013 ASPiCTaSPECL/CSJBK 42230 
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oaioi4. 

ft.} 102 ? 

001024 


C 

C 


go 1026 

00 1034 

00103/ 


ft } 1 0 4 1 

ft} 1065 

iioio/? 

■ .' Uil - 

Ooll 13_ 


J.5. .7 b 3 <» * S 2 tCL**-C SU ao I aT±» ) - 

RATlUd laSCHl T/X. 

>lATiiJ(?)='iTrtAX/A 

..ptlt.lMlNARY CJMPiJTATlOMS _ 

PMiSUal .♦ 4 . 6 *ASPLCT '*»1 » 6 b 

xxxstis ( STma x_/ jxs 1 _«y /J •J?!*' 

PHIXX = S 0 H 1 (PHISSJ) 


XKVARd )»1 ,1*SCHIT»S(5RT ( 3 . 1 41 6* 5PECU /S3RT ( ^HlSa.l , 2 *xxxS01 
XXVAR( 2 ) = 1 .l«STMAX«5r)nn3,1416»5PECL)/SaRT(»HIS3-1.2*XXXS0) 
XKVAR(3) = UI»STMAX'‘S0RT(3.U16 *SPcCU)/S3HT(ph 1S3-1.2*XXXSQ} 
.l£__y yLL ILi.Z1_gxlJJ1_215 


0)1 121 
iKl 23 _ 
001126 



DO 26' Nsl »3 
XKVAR{N)=XKVAR (N),*S0RT lo.iioll 

CHE:x ri-S STRESS_.RATio L-E 
DO 265 N=1 ,2 

IF 1 X K 'iAiU N 1_ . LT . 1 . ) R 0 TO 265 


_C 

C 


IF (RATIOINI .'iT..01 GO TO 2^0 


Oft! I3l_ 

001132 

JoiiasL 

ftnl 137 


(1.11140 



c 

c 


VALID 

XN(\|)aXN(N)*l,_ 

SDM(N)=SUM(N1 ♦xkvaR(n) 

SS U(M) =SSQ (Ml ♦XXVftR (rj> •»? 
DO TO 265 


INVALin 
_KFLAG(M)=1H* 
N0A51=i 
CQ.'ll INUL 


1T = 3 
PRI VI 
WRITE 


ONE LluE OF INPJT/OjT’UT DATA 

(6*425) lUENT«IT£'VP,TTStTJStTHICKtWipTH»C5U30»SPECL.SCRlTf ST 


001220 


001220 

041221 

61*1222 


IMAX* ( (XKVAR (N> *XFlA(>(N) ) •N sl*^) • RAT 10(1) .RATIO (?) * ASPECT* IREf 

__klfciis=LIuEid_ 

60 TO 135 


C 

27 u_ 

275 


iF\.AST (;«R0 Was blanx set ^oresZi if eof* set >iore«i 

MORE=? _____ . 

'GO to 2H0 "■ ■ 

M0R£=1 


0.t1223 


C 

C _ 
28o 


C0M3DTF AND PRINT 
'DO 265 Ns 1.3 


statistics ip Xn GT 1 


A22A0 , 
A225U 
A2260 
*22^0 
AR28() 
A2290 
_A2300 
A2310 
A2320 
A2330 
A?3A0 
A 235 y 
_A236q_ 
A2370 
_ A238u 

A2390 

A240U 
A2410 
. A2420 
A243o 

A244y 

A2450 
A246U 
A247q 
„A24Bo_ 
A2490 
A2500 
A2510 
A2520 
A?53o 
. A2540 
A2550 
A25<i>0 
A2570 

A2580 

A2590 
_ A2600 
A2610 
A2620 
A263U 
A264(< 
A265o 
A266o 
A26?0 
A2680 
A2690 


0 ;jl ?25 AVG(N)so, __ 

’001226 ' 0 £V(N)sr;, - - 

661227 IF (XN(N) .LT.2.) go to 28S 

00123? AVG(nVsSU 6 (N) /XN ('n'^ 

_6r»1234 DEv{N)sSDWT( (SSQ(N)-SU6(NV»5D'1(N)/XN(N) )/(XNl(N)-f ,) 1 

551246 285 continue 

Osll250 wRlIE (6*455) 

' 001254 IF ( (XN( 1 ) .LF.I .) .and. (XN(2> . UE. I .) .AND'.( XN (j). LE. 1 .)> '’<»() TO 31 O 

c BRA\)CH To appropriate oJTsJT format 

001273 GO TO (295.29jV36’0'W390>'39i:,396’) ♦ IT 


A?700 
A2710 
A2720 
A2730 
A274Q 
A275o 
A2760 
A2770 
A2780 
A 2 790 
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c 

':i 


001313 


wHUt (6*430) atv 

A2RI0 

6.11321 


(i(J 10 JIIJ 

A2R20 

r 6.11322 

295 

WT^nt 1694?0) AV6 ( 1 > , AVo (7) 

A2B30 

; QL>1332 ... 


wfUlE (6*435) DEv(l) .llEv (?) 

4^840 

1 0.11 3<i? 


GO ro 305 

A2950 

00]3>3 

. 30iL . 

WRITE (6*410) »VG(1)*AV5(?) 

A2B60 

0IJ13S3 


wRITF (6*44u) OEvm fDEV(?) 

A207O 

001363 

3g5__ 

.LlNis«UNEs*2 . 

A2B80 

0'ji3bS 

31U 

(i>»4S5) 

A2890 

0'ji37i . _ 


LlNES = LlNd;S^2 

A?900 

0(j1373 


(ij To (315*320) * viORF 

A2910 

P 0 1 4 u 1 

315 

IF (NOAST.E(J.l ) WRITE (6,4491 IRCT 

A2920 

onun 


IF (f'inASr.EQ.2) WRITE <6*430) 

A2930 

Q&l4l7_ 


GO ro no 

A2940 

001420 

32t> 

IF (LImE5.lt. 35) 5o Tj 124 

A2950 

Oftl423 


IF (Nn/iSl .tU.l ) WRITE (6.449) IPCT 

A?96o 

05i^32 


IF (NOiiST.EU.a) WRITE (6*490> 

A?97o 

JHoQ 

725 __ 

.lF__ .((lTYPEt5)/b,7jr, U STOP 

A2980 

goio<f 7 


NQAS f = ; 

A2990 

DU 1450 


.LIN£*5=:l* ..... . .. _ __ 

. . _ A3000 


C 

CHEC^ (.MPJT A'JO OJTPUT TY9£ ANO 8RANCH TO ARPRORITATE FORMAT 

A3010 

oSl^bO 


IF (IliNIT.EU.l) (iD To 330 

A30^0 

001^53 


IF ( 1 rt£ 7 • £Q • p ) GO TO 35o 

A3o3o 




GO TO 341. 

Alf)4o 

0 . 1 1 45s 

33u 

IF (IMFT.EQ.?) <50 TO 335 

A305o 

0^1^57 


GO TO 345 

A3o^O 

Ojl 46o 

335 

THICX=THlCX»o.o254*lnou. 

A3o70 

0 0 1 4 6 2 


GO 10 350 , 

AlOSO 

00146? 

"34 0 

THlCX = THlCx/?3**+ 

A3090 

-- 

.345 

GO ro (355»36Uf3bS*36bi385»i«5i f ITYPE 

A3100 

0.,»14 76 

350 

"(10 10 (37j*375.3df).3H0* 3B4.3H5) * ITYPE 

A3110 

OOjSlO „ 

355 _ . 

wRIlE (6.470) THICX, ITA, IT3* ITC 

A3120 

nfil 524 


Wrillr. (6*535) <Om^nT 

A3130 

00lb3? _ 


wRHE (6.500) 

A3140 

65)535 


GO 10 3B5 

A3i5o 

0J1537 

360 

wRlIE (b,.f80) ThICX,ITA,lT3,ITC 

A3160 

001553 


wHIIE (6.935) <(JME^T 

A3170 

.*> jl 561 


WRllE (6.310) ^ 

A31B0 

001565 


GO 10 3(35 

A3190 

__ 6r»1566._ 

_ 35.5 __ 

WRITE (o,.*90) THICX, ITA* IT3*ITC 

A3200 

001 60p 


WRITE (6*5j5) KOrEIMT 

A3J10 

— oi'ifeiQ 


WRITE (6*520) .. 

A322Q 

001614 


60 TO 3B5 

A3230 

0.il61S 

.37.0 

WHITE (b.tTS) ThICX*!TA*IT3*ITC 

A3240 

DDifiSi 


wHKE (6,535) KUmEMT 

*3250 

0.J1637 


WRITE (6.5C5) 

A3260 

nrii 043 


GO TO 3B5 

A327o 

d51‘>44 

. 375 .. 

WRITE (b.tHS) THICX, ITA* IT3* I TC 

A3?8o 

0 1 0 6 0 


wRiic (6*535) AOm-NT 

A3290 

O.j 1 666 


write (6,515) 

A3300 

"6.0167? 


GO To 3B5 

A3310 

651673 

38V„. 

WRITE (6*495) TH 1 CX * 1 T A * 1 73 * I TC 

*33^0 

Oil 1 707 


WRITE (6*535) KOMiiMT 

*3330 

6r»171S 


WRITE (6*525) 

*3340 

6.11721 

38i> 

IF (M0RE.EQ.2) GO TO 125 

*3350 
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t.c. 


fi A i 7 ? 



-0-0 1 / 

"531 7c4 
ftOl 730 

c 

c 

tRRDW MESSAGE - MTSSIMG DATA 

HKIMl 5Su 

GO 10 135 - .. 

oTiuai 

395 

C 

CONTINUE 

c 

c .. ^ 

0M731 

40U 

format (lX,A8»l‘ff^X,F5,lf3XfF3,if2XtK6.i,3AfFb,3,3A»F&,3f3AfFb,3»2 
lX»Fb*3''^^*^5.2fAi*F8,2’^l»*b.l»fb.liFBt3»3XfF5*3»A7) 

Alii 731 


format (jX,A«,l4,iX,pX*X5.j,jX.rb.i«iX*Fb.3*2X»Fi>.,2,2X»F6,2*2XiF6» 

12»3X,Fb.2,2XjFb.2-» 2X iFb.,2 , A 1 , 1 X ? K6. 2 , Al , J X , F6 »AlAJ.iLX t.F6,,3j 2Xlf 


nf*l73i 

DiJl’ii 

9..1731 

0^1731 


. <►» W- 

4i3 

>2\». 

42& 


nnl73i 

001731 

P(ji731 
_00i73u 


43u 

>35 

44t> 

_44S> 


6.3i A7) 

_ F 0 ><<tAT . (64 X,ohVERA 3E VftLUr«»-9,2»5XiF7.2) . . . 

FOH'IAT (6)Xi«AvE,kA3F vAi,U"<*i3X*Fo«2»2’^*^^*2*2**^*^*2* 

.. Fy3v|AT . (59X ,ttAVtMAt.t; V nUU t «* » 4X, Fo, 2 , 3X * F(,, 2) 

FUR'^at nx,AR,lX,I4,F8.1,r6tI,Fd.3,FH.2,F8.2,F«J»?,F8»2,F8.2,FH«2,A 

.-1JiA.xjs.fi, 2 s ai , xy . 3 • F S ,3 s ^8. 3 » 2 X ♦ a 7| 

FORMAT (b1 X,»sTAN)ARD U£ V . » , 3 X , Fo .2 » 2 X » F (>.2 t 2 X ,Ffe. 2 X ) 

_ Format (b'^X»* 3 lAM 3 AkO 0EV.»iAXtFb,2»3x.F6.2) „ . 

FORMAT 164X,«STAn3AH0 OE V . «• 9. 2 ,5X, F7 ,2) 

format ( 2 Ho«,'*MOTE- net S-criJ\| stress exceeds *,13** percent 


A33&U - 

A3370 
A3380 
A339J 
A3A00 
A3AlO 
A3A20 
A343U 
A3A4o 
A3A5o 
A3480 
A 34 ^0 
>3480 . 
A349U 
A3500 
*3510 
A3520 
A3530 
A35A0. 
A3580 
A35(>0 . 
A3570 
A 35 OO 
A3590 


-8U1711 A5Jl 


OF' 

iiEO strength, value Not inclujeo in avg, value, or std. dev. ♦) 

format (2H n«.ttN0TE" TNVALT3 MC. yALU£i_NOI,.lNCi.JuiO_IN AVGt— VALUE A3600. 

A36IO 


001731 

91*11731 

-OOlTBl 

0>*173i 


ASS. 
46u 
A65 , 
470 


lOR 5T0. 0£V.») 

Format (2A7) . , - 

format (*PM if 0I5P0SE, PRINT JNLINEU ON 512. PH«» 

FORMAT (3X,Ii,/X,AB,2A13,5X,2»XiF5.3) 

format (♦)Tab,e <c and kjC fracturei toughness Values for 

ICOnRAC T TY PE SPECIimEnS QF »,F5«3»* ImCH THICK _*A9, AlOi* AllOY . * t AJ- 


oji73i 

-^75 

p(;i731 

1 

_ ._48 a 

niTl73i 

485 

nai73L 

4_9u 

i.rvii?jL 

4?5.__ 

001731 

500 


20 /) 

format (*iTAR_E 

iSPLCIMENS OF tt,F5. 

Format (oitar.e 
luGHNEsS OF *,rS«3tA 


FRftCIJRE toughness VALUES ?DR COMPACT TYPE 
2»* MIL.iIMETER THICK *A9iAiq*» ALLOY *Aiq/) 


PJANE-STRESS ANJ 
INCH thick **a 9 »A 1 o»» 


TRANSITtDNAL fracture to 

ALkOY *,AlO/) 


-flO 17-31-.- 5a5. 


FORMAT (MtTAH.E PLANr-STRESS AND TRANSIT ION«.t FRACTURE TOUGHN 

7 eVs 0 F»^r 5 '. 2 * A MILLT^ £R hick *A 8 » Aiu»*ALlOV»A, 6 /> 

FoRm*T (mITAH.E ^ SjR"ACE-FLAW fracture TouSHnESS oF *,F 5 . 3 ».. 
i« INCH THICK »,A 8 *Alfl»* AiJlOY »»A 1 ()/) 
format t*lTAR_E SURF aCEhF|_aw FRACTURE T 0 U 0 Hn£ 6 S OF *,F 5 . 2 ,* M 

ULLIiMETER THTCK *,A 8 , AlOtttA.LjY **Al 0 /) 

FORMAT {IX,*SRECIMEJN TEST _m** 1 ENS 1 LE**+ ♦♦♦♦Sp-CIMEN D 1 ME,nSI 0 NS_. 

^;,j^HNESS SI 2 Ei RE 3 T LOAD STRESS 

TE'> RRUPERTlcS THICK- wIOTh CRACK 

OFFSrT APPRNT OFF- APP- RATIO RATIO* 

length * 39 X.*SET RNT »/ 12 X,»F. TYS TO 

{A)» 23 Xt«K(IO K(APP) *» 8 X»* PQ/PmAX S 

INC-Hi I.nCH KIpS 

*//> 


!♦♦♦ ♦♦♦♦lOa3**m* 

2 _REF*/2X,* IUEnT 

3 P(Q) P(MAX) 
4/34X»*fjESS 


55 (1) (w) 

bN /iYs */i7x.»Ksi <si Inch , 

/kKPS ♦♦kSI-SURT t inch) ♦♦ INCH INCH 
format (1 X.*S=’ECIM£N test ♦♦♦TtNSILE**M ♦♦♦♦50-CIMEN DIMENSIONS 
(♦♦♦ ♦♦♦♦LOa 3«*** MAfERlA-l toughness SI2£; REUT LOAD STRESS 

2 REF*/2x,* lOENT TEM» PROPERTIES IMICK- iHlOTH CRACK 

3 R(i)) P(MAX) OFFSrT APPRNT OFF” APP" RATIO RATIO* 

4/34x,*ness length *n9x,*set rnt » /i2X, * c ns TU 

5S (f)' (W) ■ (4ri22Xf*K(IC) K(App) 


*5X,*PQ/p 


A3620 

A3630 

A 36 A 0 

A 365 O 

A3600,- 

A3670 

A360u 

A369o 

A370O 

A37lO 

A3720 

A373O 

A37A0 

A37SU 

A3760 

A3770 

A37B0 

A37VO 

A3800 

A3BI0 

A3320 

A3830 

A3S40 

A3850 

A38bO 

A387C 
A3B80 
A389U 
A39OO 
' A3910 
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■) 




001731 513 


— SU/IYS s /1 S* + M\|/3J (M) ♦♦♦ ♦♦♦♦♦♦♦♦'1lLLI'1ETER***** + t K 

TlLO'iE.VTUNS** ♦♦♦♦MN/ ('ll MILLIMETER *//) 

FijRvlAT (M,»s = KC1mEh test TEnSIlE ♦♦♦♦♦SrEcIsEn UImEnSIONS** 

l*t*» ♦♦♦♦STP-5S+ + ** ♦♦^AIEPlftL TOObH'lESS^* ♦♦♦NET STrtEisS RATIO 

, 2** REF */ 2 x»*ID-NT PR0PERIIE5 TMlCK- wlOTH ♦ C«aCK LE 

jNtiTi ♦ OFFSET MAXI- OFFSET ARPHNT CHITCL S ( -) /TTS 1 1 -?C « • ) / 

j4wiimj<a±_yJE^iL.Utl N&S5 InUIAj E!1ni«l. - - - mU 

OFFSET APPRNT CR1TCL*/12X,*F STR 

.6 STR._ . (T) (W>_. ?CO> 2C(C) S(0) S(C» K(Q1 

7K(A3R) Mc)*/i2X*tt <51 KSI INCH INCH INCH INC 

_8H KSI KSl ♦♦♦♦<SI-50RT(I.NCH) ♦♦♦ »//> 

format (1X***S3£CImEN tEst tensile ♦♦♦♦♦SrEcIMEN 0IMEN5T0Ns^+ 

!♦*♦* ♦ ♦♦♦SIRESS^^^^ ♦ ♦material TOUOHNESS^^ ♦♦♦NET STRESS RATIO 

REF »/2Xt*lDENT TfM^i PROPERTIES HICK- WlUTH ♦ CHaCK' LE 

JNGH ♦_ OFFSET MAXT^ OFFSET APPRNT CRUCL S (•) /TYS I1-2C (-> / 

A«)*/16X,* YIFLO ULT NEsS INITIALI FuNAL MU 

.5*’ . ... ... OFFSET APPRNT CRITCL*/11X,* C SR 

6 STR (T) (W) ?C(0) 2C(C) SO) SIC) 


— oiiiJl] 52 


Ji»I73I 


.0'n73l_ 

0*0731 

.5i')T3l„ 
010 731 
_001731 
0.0 731 


, - K(Q) 

7 K ( App) ) K(C)*/1_2X«0 MNI/SjIM) ♦♦♦♦♦♦♦♦♦♦mILLIMETEr^^^^^^^^ 

({♦♦ ♦♦♦mn7sQ(M) ♦♦♦ ♦♦♦♦♦♦♦MN/M3/2*+^*^*^*/7) 

FJRMAJ (|X,»s3£CIMEn test ♦ ♦ ♦ I EnSILE^ ♦♦ ♦♦♦♦♦♦SPECIMEN DImENSIO 
IfjS^^*** ♦♦♦♦STRESS^^^^ material TOUGHNESS ♦♦NET STRESS^^ A/ 

22C rEF*/2X,*I0EmT TEsI PRqpErTIES THICk- WIDTH ♦♦CrACK 
3 OFFSET MflXi- offset A^PRNT RAT I0»/16Xf Yl 

AELO t.|LT NESS LEMUTH DEPTH _ MDM* 26 X**PFF 

APPRNT**/12X,*F STR STw" (f) (In " 2(C) A( 

S(0) SIC) „ K(0) K (APP) o/ieXfAKSi KSI INCH 

INCH TNCH inch KSI KSI ♦♦KSI-SORT ( INCH) ♦♦♦♦// ) 

-525 Format ( 1X,*s3ECTMEw jEsr ♦♦♦TENsILE^^^ ♦♦♦♦♦♦specimen UImENSIq 

♦♦♦♦STRESS^^^+ MATERIAL TOUGHNESS ♦♦NET SIRESS^^ A/ 

22 C. _..WFF«/ 2 X,<*IUENJ lEsI ^RO.PEHJI£S_._..THJCkl-.:_ width ♦♦CRACK 

3 SI/E^** offset MaXI- offset A^PRNT RATIOA/lbXfO YI 

AELI) ULT :^^ESS LEnGTH OEpTH__ _ MUM*26xt*0FE 

5SET APPRNr»/i 2 X,ttC SIR STR (T) (v) 2(c) A( 

60) S(0) 5(C) K(Q) K(APP)»/l7X»*^^^viN/S3(M)^^ ♦♦♦♦♦ 

^♦♦♦♦♦•MII.L1HETER^^^^^^^^^^ ♦♦♦MN/Sa(M) ♦♦ ♦♦♦♦MN/ (M) o/p^^^^t*//) 

F0RMAT_J6A10) 

(a5x,baio/7) 


5SEI 

60 ) 

7 


53u 


535 

5 ao 

5hS 

.5.5y_ 


format 

FORMAT 

FoR-mAT 

format 

END 


(HX*FS«U»5 ;'X»ti) 

(Aa,u"F5.u. lI.A2,3x.U,A7) 
.<*. *“-".J.M.1.S_LI!NE OMITTED - 


MISSING T YS OR JPEC TYPE*) 


A392U. 
A 393 O 
A39M0 
A3950 
A3960 
A 397 O 
_ A3980 
A3990 
AAOOO 
AAolO 
AA020 
A4030 

aaoao. 

aaoSq' 

A4060 

AA07u 

AAQUO 

Aa09U 

A4I0U 

AAllo" 

AAi 2 o 

Aa 13U 

AAl AO 

AAISU 

A 4 I 6 U 

AA170 

AAiao 

AaHU 

AA200 

AA210 

Aa 220 

Aa 2 J 0 

A42A0 

AA250 

AA 2 G 0 

AA270 

AA200_ 

AA290 

AA 3 OO 

AA310 

AA320 

AA33o« 


173 


APPENDIX H 


THE WEIGHTING FUNCTION, W(X) 

In order to obtain uniformity of variance in combined sets of fatigue data, a 
weighting function, W(X), was applied to each data population. The following 
comments are included to clarify the method employed in this weighting process. 

Initially, the data were analyzed without weights, and the quality of fit 
was based on the parameter. The process used in maximizing involved a 
minimization of the sum of squares of deviations or SSD. Using equations (2) 
and (4), the unweighted SSD can be described as 

n _ n _ 

SSD = 2 (Y. - Y)2 = 2 (Y. - Ao - A,X)^ (HI) 

i=l ^ i=l ^ 

where X represents the mean, or predicted value of X for a particular value of Y . 

1 

When the weighting function was used, the minimization was based on a modified 
SSD, written as follows, 

n 

SSD = 2 W. (Y. - A^ - A,X)2 . (H2) 

i=l ^ 

In this way, each deviation from the mean is modified according to the magnitude 
of W^. When W^ is small, the square of the residual is reduced correspondingly. 

If W^ is near unity, almost no modification of the residual results. 
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